EXPERIMENTS WITH ZETA ZEROS AND PERRON'S FORMULA 



ROBERT BAILLIE 

Abstract. Of what use are the zeros of the Riemann zeta function? We can use sums involving 
zeta zeros to count the primes up to x. Perron's formula leads to sums over zeta zeros that 
can count the squarefree integers up to x, or tally Euler's (j) function and other arithmetical 
^ functions. This is largely a presentation of experimental results. 
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1. Introduction 

Most mathematicians know that the zeros of the Riemann zeta function have something to 
do with the distribution of primes. What is less well known that we can use the values of the 
zeta zeros to calculate, quite accurately in some cases, the values of interesting arithmetical 
functions. For example, consider Euler's totient function, (f){n), which counts the number of 
integers from 1 through n that are relatively prime to n. The summatory function of (p, that 
is, the partial sum of (f){n) for 1 < n < x, 

n<x 

is an irregular step function that increases roughly as x^. In fact, a standard result in analytic 
number theory is that [13, Theorem 330, p. 268] 

Sx"^ 

^(x) = — 5" + 0(x logx). 

But such an estimate mainly tells us how fast increases. ^ is a smoothly increasing 

function and does not replicate the details of the step function at all. But by using a 
sum whose terms include zeta zeros, we can get a function that appears to closely match the 
step function ^{x). Our approximation rises rapidly at the integers, where has jump 

discontinuities, and is fairly level between integers, where ^{x) is constant. 

As another example, consider a{n), the sum of the (positive) divisors of n. The summatory 
function 

T{x) = J2^{n'). 

n<x 

increases roughly as x^, but zeta zeros give us a good approximation to T{x). Both of these 
sums come to us courtesy of Perron's formula, a remarkable theorem that takes integrals of zeta 
in the complex plane and yields approximations to the sums of integer arithmetical functions. 
We then use Cauchy's theorem to estimate these integrals by summing residues; these residues 
are expressions that involve zeta zeros. 
Asymptotic estimates of the type 

T{x) = fix) + 0{g{x)) (1.1) 

are known for some of the summatory functions we consider here. See, for example, section 25. 
It often turns out that the method used here (summing residues) easily produces 

T{x) ~ f{x) + sum over zeta zeros, (1-2) 

and even gives us the exact coefficients in f{x), although it does not yield the error estimate 
g{x). Assuming that equation (1.2) holds, equation (1.1) is essentially estimating the size of a 
sum over zeta zeros. 

Using zeta zeros, we can even do a pretty good job of computing the value of vr(x), the number 
of primes < x, although Perron's formula is not used in this case. But all these examples suggest 
that zeta zeros have the power to predict the behavior of many arithmetical functions. 
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This paper is long on examples and short on proofs. The approach here is experimental. 
We will use standard results of analytic number theory, with appropriate references to those 
works, instead of repeating the proofs here. In many cases, we will give approximations to 
summatory functions that are not proven. Proofs that these expressions do, in fact, approach 
the summatory functions could make interesting papers. Further, this paper is still incomplete, 
and should be considered a work in progress. 

Many thanks are due to professors Harold Diamond and Robert Vaughan, who have patiently 
answered numerous questions from the author. 

2. The Riemann Zeta Function and its Zeros 

Here is a brief summary of some basic facts about the Riemann zeta function and what is 
known about its zeros. As is customary in analytic number theory, we let s = a + it denote a 
complex variable. 3fJ(s) and are the real and imaginary parts of s. 

The Riemann zeta function is defined for 3fJ(s) > 1 by 

oo _ 
n=l 

Zeta can be analytically continued to the entire complex plane via the functional equation 
[34, p. 13], 

C(.) = 2%-iC(i-^)sin(y)r(i-.). 

Zeta has a simple pole (that is, a pole of order 1) at s = 1, where the residue is 1. 

The real zeros of zeta are s = —2, —4, —6, .... They are often referred to as the "trivial" 
zeros. There are also complex zeros (the so-called "nontrivial" zeros). These zeros occur in 
conjugate pairs, so if a + bi is a zero, then so is a — 6i. 

When we count the complex zeros, we number them starting at the real axis, going up. So, 
the "first" three zeros are approximately 1/2 + 14.135i, 1/2 + 20.022i, and 1/2 + 25.011i The 
imaginary parts listed here are approximate, but the real parts are exactly 1/2. Gourdon [10] 
has shown that the first 10^^ complex zeros all have real part equal to 1/2. The Riemann 
Hypothesis is the claim that all complex zeros lie on the so-called "critical" line Re{s) = 1/2. 
This is still unproven. However, we do know [12] that the zeta function has infinitely many 
zeros on the critical line. We even know [6] that over 40% of all of the complex zeros lie on this 
line. 

3. Zeta Zeros and the Distribution of Primes 

Let c be the least upper bound of the real part of the complex zeta zeros. We know that 
1/2 < c < 1. Then it follows that [21, p. 430] 



7r(x) = li(x) + O (x''log(x)) 
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where 




1 



dt. 



2 



log(t) 



This shows a connection between the real parts of the zeta zeros and the size of the error 
term in the estimate for '7r(x). However, this sort of estimate says nothing about how 7c{x) is 
connected to the imaginary parts of the zeta zeros. We will see that the imaginary parts of the 
zeros can be used to calculate the values of 7r(x) and many other number-theoretic functions. 

This topic falls under the category of "explicit formulas" . Two examples of explicit formulas 
are commonly found in number theory books: a formula for 7r(x), and a formula for Chebyshev's 
psi function. We will examine these, then consider many more examples. 



von Mangoldt expanded on the work of Riemann and derived a remarkable explicit formula 
for 7r(x) that is based on sums over the complex zeta zeros. 

First, note that, if x is prime, then n{x) jumps up by 1, causing 7r(x) to be discontinuous. 
7r(x) is continuous for all other x. 

For reasons that will become clear in a moment, we will define a function vro(x), to be a slight 
variation of tt{x). These two functions agree except where 7r{x) is discontinuous: 



When 7r(x) jumps up by 1, 7ro(x) equals the average of the values of 7r(x) before and after the 
jump. We can also write vro(x) using limits: 



4. The Riemann-von Mangoldt Formula for 7r{x) 




n{x) — 1/2 if x is prime, 
n(x) otherwise. 




The Riemann-von Mangoldt formula [27, pp. 44-55] and [28, eq. 5] states that 




(4.1) 



n=l 



where 




rft-log(2). 



(4.2) 



Here, is the Mobius mu function. If the prime factorization of n is 



then /i(n) is defined as follows: 
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1 



if n 



1 







if any ai, 02, ■ ■ ■ , 
otherwise. 



ttk > 1, 



(4.3) 



(-1) 



That is, /i(n) = if n is divisible by the square of a prime; otherwise is either —1 or +1, 
according to whether n has an odd or even number of distinct prime factors. 

Several comments about the Riemann-von Mangoldt formula are in order. 

First, notice that the sum in (4.1) converges to tto{x), not to 7r(x). That's why we needed to 
define 7ro(a;). 

Second, the sum in (4.1) is not an infinite sum: it terminates as soon as x^/" < 2, because 
as soon as n is large enough to make x^^" < 2, we will have /(x^/") = 0. 

Third, the summation over the complex zeta zeros is not absolutely convergent, so the value 
of the sum may depend on the order in which we add the terms. So, by convention, the sum is 
computed as follows. For a given A^, we add the terms for the first N conjugate pairs of zeta 
zeros, where the zeros are taken in increasing order of their imaginary parts. Then, we let N 
approach 00: 



where the asterisk denotes the complex conjugate. Here is an equivalent way to write this sum: 



where t is the imaginary part of the complex zero p. 

Fourth, as a practical matter, when we compute the sum over zeta zeros, we need to use only 
those zeta zeros that have positive imaginary part. Suppose that a+bi is a zeta zero with 6 > 0, 
and that the corresponding term in the series is A + Bi. When evaluated at the conjugate zeta 
zero, a — bi, the term will have the value A — Bi. These terms are conjugates of each other, so 
when they are added, the imaginary parts cancel, and the sum is just 2A. Therefore, we can 
efficiently compute the sum as 



where now the sum is over only those zeta zeros that have positive imaginary part. 

Can (4.1) and (4.2) really be used to approximate the graph of 7r{x)7 Yes! Figure 1 shows 
the step function 7r{x) and some approximations to tt{x). To make the left-hand graph, we 
used 10 pairs of zeta zeros in the sum in equation (4.2). The step function 7i{x) is displayed as 
a series of vertical and horizontal line segments. Of course, the vertical segments aren't really 
part of the graph. However, they are useful because they help show that, when x is prime, (4.1) 
is close to the midpoints of the jumps of the step function. 





00 



2^^(5^h(x'"=)) 



k=l 
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In the right-hand graph, we used the first 50 pairs of zeta zeros. For larger values of x, this 
^ives a better approximation. 



15 - 



10 




15 - 



10 




Figure 1. approximating 7r(a;) with 10 and 50 pairs of zeta zeros 

In Figure 2, we have plotted (4.1) for 100 < x < 130, using 50 and 200 pairs of zeta zeros 
in the sum in (4.2). This is an interesting interval because, after the prime at x = 113, there 
is a relatively large gap of 14, with the next prime at x = 127. Notice how the graphs tend to 
level off after x = 113 and remain fairly constant until the next prime at x = 127. It is quite 
amazing that the sums (4.1) and (4.2) involving zeta zeros somehow 'know' that there are no 
primes between 113 and 127! This detailed tracking of the step function 7r(x) (or 7ro(x)) is due 
solely to the sums over zeta zeros in (4.2). 




120 125 



130 



Figure 2. using 50 and 200 pairs of zeros to approximate 7r(x) for 100 < x < 130 



Stan Wagon's book [36, pp. 512-519] also has an interesting discussion of this topic. 
In section 27, we see how to use zeros of Dirichlet L-functions to count primes in arithmetic 
progressions. 
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Figure 3. using 1000 and 2000 pairs of zeros to approximate n{x) near the jump 
at X = 17 




Figure 4. Partial sums of 10 and 40 terms of a Fourier series with unit jump 



5. The Gibbs Phenomenon 

When a Fourier series represents a step function, the ra*^ partial sum of the Fourier series will 
undershoot and overshoot before and after the jumps of the step function. If the step function 
jumps up or down by 1 unit, there will be an undershoot and an overshoot in the infinite sum, 
each in the amount of 

1 1 r sin(t) , , , 
g = — + - / — — dt ~ 0.0894898722 (5.1) 

2 TT Jo t 

before after the jump. This is called the Gibbs phenomenon. Figure 4 shows the overshoot in the 
partial sum of the first n terms of a Fourier series that has unit jumps at a; = 7r/2, x = 37r/2, 
X = 57r/2, etc. In this figure, the first peak after the jump occurs at about x = 1.72788, 
y = 1.08991 for n = 10, and at about x = 1.61007, y = 1.08952 for n = 40. (For any given n, 
the first peak occurs at x = | + ^). If the step function has a jump of c, then the undershoot 
and overshoot are both c times the constant in equation (5.1). 

Figure 3 shows a closeup of the approximation (4.1) to 7r(x) near x = 17. (There's nothing 
special about 17). This looks very much like the partial sum of an Fourier series for a step 
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function, including a Gibbs-like phenomenon before and after the jump. In the right-hand side 
of Figure 3, which is a sum over 2000 pairs of zeta zeros, the minimum value just before the 
jump is about 6 — 0.08885; the maximum value just after the jump is about 7 + .08983. This 
overshoot is visible in the some of the graphs in the sections below. 

6. Rules for Residues 

Later, we will need to compute the residues of functions that have poles of orders 1 through 
4. Formulas for poles of orders 1 through 3 can be found in advanced calculus books (see, 
for example, Kaplan [16, p. 564] for the formulas for orders 1 and 2 and Kaplan [16, p. 575, 
exercise 8] for the formula for order 3). For convenience, we state these formulas here. Their 
complexity increases rapidly as the order of the pole increases. However, Kaplan [16, p. 564] 
also gives an algorithm for computing the residue at a pole of any order. Instead of writing out 
the formula for order 4, we will give a general Mathematica procedure that implements that 
algorithm. 

If f(z) = where A{z) and B {z) are analytic in a neighborhood of s, and where A{s) ^ 
and B{z) has a zero of multiplicity 1 at z = s, then the residue of f{z) at z = s is 

re3(/W,3) = ;^. (6.1) 
If the conditions are the same as above, but B{z) has a zero of multiplicity 2, then the residue 

is 

, , QA'(s)B"(s) -2A(s)B^^\s) , , 

If the conditions are the same as above, but B{z) has a zero of multiplicity 3, then the residue 

is 

120B^^\sfA"{s) - 605(4)(s)5(3)(s)A'(s) - 12A{s)B'^^\s)B^^\s) + 15A{s)B^^\sf 

405(3) (s)3 ■ ^^-^^ 

Here is Mathematica code that computes the formula for any order m, m > 0. 

residueFormula[m_Integer?Positive, z_Symbol] := 
Module [ 

(* compute the residue of A[z]/B[z] where 

B has a zero of order m and A != *) 
{ aa, a, bb, b, k, ser, res }, 

aa[z] = Sum[a[k] z"k, {k, 0, m - 1}]; (* need z"m through z"(m-l) *) 
bb[z] = Sum[b[k] z"k, {k, m, 2*m - 1>] ; (* need z'm through z"(2*m-l) *) 
ser = Series [aa [z] /bb [z] , {z, 0, -1}]; (* need z'm through z"(-l) *) 
res = Coefficient [ser , z, -1]; (* residue = coefficient of z"(-l) *) 
(* replace a[ ] and b[ ] with derivatives *) 

res = res /. {a[k_] -> D[A[z] , {z, k}]/k!, b[k_] -> D[B[z], {z, k>]/k!}; 
Simplify [res] 

] 

For example, the following command will generate the expression in residue rule (6.2): 
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residueFormula[2, s] . 

In various sections below, we will first define A{s) and B{s), then use residueFormula to 
expand the result. We do this, for example, in section 17. 

Finally, later (sections 13 and 17), we will need a formula for the residue at z = 1 of 
expressions like 

CM! 

C(2z)- 

where the residue (here, of order 2) arises not from a zero in the denominator, but from a pole 
in the numerator ai z = 1. If f{z) has a pole of order n at z = s, then the residue at z = s is 
[16, p. 564] 

J-n— 1 ^ ' 



(n — 1)! z-^s dz 

7. Perron's Formula 
Given a sequence of numbers a„ for n = 1, 2, . . . , a Dirichlet series is a series of the form 



E 



On 



n=l 

where s is a complex variable. If a Dirichlet series converges, then it converges in some half- 
plane Re{s) > cq. 



Now, suppose we have a Dirichlet series 



On 



F{s) (7.1) 



n=l 

that holds for Re{s) > cq. 

Perron's formula inverts (7.1) to give a formula for the summatory function of the 
function of x: 



E 



1 



c+iT 



an= lim — / F{s)—ds. (7.2) 

n=l ^ c 11 

This holds for x > and for c > cq. Whenever we see an x as the upper limit on a summation, 
we mean that the sum extends up to [x\ . 

The prime on the summation means that, when x is an integer, instead of adding to the 
sum, we add ^a^. Let G{x) be the sum on the left-hand side of (7.2). Suppose x is an integer 
and that > 0. Then G{x — .1) does not include the term ax at all, G{x) has as its final 
term, and G{x + .1) has ax as its final term. Thus, at x, G has a jump discontinuity, and the 
value of G{x) is the average of, for example, G{x — .1) and G{x+ .1). Of course, this argument 
also holds with .1 replaced by any e between and 1. 

To summarize, the left side of (7.2) is a step function that may jump when x is an integer. 
If this step function jumps from a to a new value b, the right side converges to the midpoint, 
(a + b)/2, as T approaches infinity. The important point is that the integral on the right side 
of (7.2) equals the sum on the left, even where the sum is discontinuous. 
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We will estimate the integral in (7.2) by integrating around a rectangle in the complex plane. 
Recall that, under appropriate conditions, the integral around a closed path is equal to 2'7ii 
times the sum of the residues at the poles that lie inside that path. We will use the sum of the 
residues to approximate the line integral. 

In all of our examples, the will be integers. It is quite remarkable that the line integral 
(7.2) in the complex plane can be used to sum the values on the left side of (7.2). It is not at 
all obvious that a complex line integral has anything whatever to do with sums of the integers 
a„. Nevertheless, Perron's formula is proved in [21, pp. 138-140], among other places. For the 
reader who is still skeptical, we will give numerical examples illustrating that the left and right 
sides of (7.2) can, indeed, be quite close when a modest number of poles is taken into account. 

Perron's formula was first published by Oskar Perron in 1908 [23]. 

The plan for most of the rest of this paper is as follows. In each section, we will apply 
Perron's formula to an arithmetic function whose Dirichlet series yields a sum over zeta zeros. 
This means that the Dirichlet series must have the form 



where G{s) has factors of the zeta function. The point in having ({s) (or C{2s), etc.) in the 
denominator is that, each zeta zero will give rise to a pole in the right-hand side. When we 
sum the residues at these poles, we get one term in the sum for each zeta zero. 

There are two sources we will frequently cite that conveniently gather together many such 
arithmetic functions: [9] and [20]. The former lists the Dirichlet series in a convenient form. The 
latter derives them, along with many asymptotic formulas. Below, we will work with some of the 
most common arithmetic functions. These two references also discuss many more arithmetic 
functions which we do not consider, such as generalizations of von Mango Idt's function and 
generalizations of divisor functions. 

In some cases. Perron's formula leads to a sum of terms involving zeta zeros, which, like 
Perron's integral, is provably close to the summatory function. In many cases, the sum appears 
to be close to the summatory function, but no proof has been written down. Finally, in other 
cases, the sum over zeta zeros appears to diverge. 



We will now present one of the best-known applications of Perron's formula. This example 
is discussed in detail by Davenport [7, pp. 104-110], Ivic [15, p. 300], and Montgomery and 
Vaughan [21, pp. 400-401]. 

Chebychev's psi (ip) function is important in number theory. Many theorems about iplx) 
have a corresponding theorem about 7r(x), but ip is easier to work with. ip is most easily 
defined in terms of the von Mangoldt Lambda function, which is defined as follows: 




8. Using Perron's Formula to Compute tp{x) 




power of a prime (where k > 0), 
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V^(x) = ^A(n). 



n=l 



It is to be understood that the upper hmit on the summation is really the greatest integer 
not exceeding x, that is, [a;J . 

One can easily check that ip{x) is the log of the least common multiple of the integers from 
1 through X. For example, for n = 1, 2, . . . , 10, the values of A(n) are 0, log(2), log(3), log(2), 
log(5), 0, log(7), log(2), log(3), and 0, respectively. The sum of these values is 

^p{10) = 3 ■ log(2) + 2 ■ log(3) + log(5) + log(7) = log(23 . 3^ ■ 5 ■ 7) = log(2520) 

where 2520 is the least common multiple of the integers from 1 through 10. 



50 h 



40 



30 - 



20 



10 



Figure 5. The graph of iIj{x) 
Figure 5 shows a graph of ip^x). 

The Dirichlet series for A is [13, Theorem 294, p. 253], [7, p. 105] 



E 



n=l 



A(n) _ Cis) 



(8.1) 



This series converges for 3ft(s) > 1. 

When we apply Perron's formula to this Dirichlet series, we will get a formula that approxi- 
mates the corresponding summatory function, iIj{x). But recall that, in Perron's formula (7.2), 
the summation includes half of the last term if x is an integer that causes the step function to 
jump. So, Perron's formula actually gives a formula for ipo, a slight variant of ip: 




if X is the k^^ power of a prime (where k > 0), 
otherwise. 



That is, where ipi^) jumps up by the amount A(x), i^o^x) equals the average of the two values 
ip{x) and ipi^ ~ !)• Recall that we made a similar adjustment to 7r(x) and called it 7ro(x). 
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Given the Dirichlet series (8.1), Perron's formula states that 

1 C'(s) 

Mx)= lim / -^—ds. (8.2) 

The computation of ipoi^) (o^' '^{^)) iiow reduces to computing a hne integral on the segment 
from c — iT to c + iT in the complex plane. Here is our plan for computing that integral. First, 
we will choose c. Then, we will form a rectangle whose right side is the line segment from 
c — iT to c + iT. Next, we will adjust the rectangle to make it enclose some of the poles of our 
integrand. (We must make sure that the rectangle does not pass through any zeta zeros.) The 
left side of the rectangle will extend from, say, a + iT to a — iT. We will integrate around this 
rectangle in the counterclockwise direction. 

The integral around the rectangle is 

pc+iT i-a+iT i-a-iT ^ fC-iT 



I{s)ds + — I Iis)ds + — I Hs)ds+—I I{s)ds. (8.3) 
where the integrand I{s) is 



The exact value of this integral can be computed easily: it is the sum of the residues at the 
poles that lie inside the rectangle. 

It will turn out that the last three integrals in (8.3) are small, so we have, approximately, 

-iT 



iI)q{x) ~ / I{s) ds ~ sum of residues. 

2vr2 J^_iT 



Our task now is to locate the poles of the integrand, determine the residues at those poles, 
then adjust the rectangle so it encloses those poles. 
First, the integrand has a pole at s = 0. 
Second, it is known [13, Theorems 281 and 283, p. 247] that 



and that 



CW = ^ + o(i) 



C'W = ^ + o(i), 



So 



C'(^) 



+ 0(1). 



as) s 

Therefore, the integrand has a pole at s = 1. 

Finally, the denominator is at every value of s such that ({s) = 0, that is, at each zero of 
the zeta function. We will see that every zeta zero, both real and complex, gives rise to a term 
in a sum for ip{x). 

Next, we list the residues at these poles. At s = 0, the residue is — = — log(27r) ^ 
-1.837877. 
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At s = 1, the residue is x. 

To compute the residue of this integrand at a typical zeta zero, we will use residue formula 
(6.1). But to use this formula, we must assume that if p is a zero of zeta, then it is a simple 
zero, that is, a zero of multiplicity 1. All of the trillions of complex zeta zeros that have been 
computed are, indeed, simple. So are the real zeros. In any case, since our goal is to get a 
formula that involves a sum over only a modest number of zeros, this assumption is certainly 
valid for the zeros we will be considering. 

To get the residue of (8.4) at s = p, we let B{s) = ({s) and let A{s) equal everything else in 
the integrand (8.4). Thus, A{s) = —('{s)x'^/s. With these choices, A{s)/B{s) is equal to the 
integrand. The residue at the zeta zero s = p will be 

Mp) _ Cip)x'' 1 _ xP 
nP)~ ^W)~~~P' 

Notice that the above argument holds equally well for both complex and real zeta zeros. 

We will now adjust our rectangle to make it enclose some of the above poles. 

First, we'll make sure that the rectangle around which we integrate encloses the poles at 
s = and s = 1, where the residues are and x, respectively. We can do this, for example, 

by moving the right boundary of the rectangle to the line 3^^(s) = 2. 

Next, we'll adjust the left-hand boundary of the rectangle so it encloses the first M real zeros 
—2k = —2, —4, ■ ■ ■ , — 2M. The sum of the residues at the corresponding poles is 



E 



-2k 

k=l 

We also adjust the top and bottom boundaries of the rectangle to make it enclose the first 
pairs of complex zeros. Note that, for a given A^, there will be a range of values of T such that 
there are A^ zeros whose imaginary parts, in absolute value, are less than T. We can choose A", 
or we can choose T, whichever is more convenient. 

The k^^ pair of conjugate zeros gives rise to the two terms 

Pk 

(where the bar denotes the complex conjugate), and 



pk 

It is easy to verify that these terms are conjugates of each other. When we add these terms, 
their imaginary parts will cancel and their real parts will add. Therefore, the sum of these two 
conjugate terms will be 
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The sum of the residues at the complex zeros of zeta will be 

k=i 

When we take into account all of these poles, the sum of the above residues is 

We can extend the left-hand side of the rectangle as far to the left as we wish, enclosing an 
arbitrarily large number of real zeta zeros. The second sum converges quickly to ^ log(l — ^) 
as M approaches oo. It can be proven that it is also valid to let approach oo in (8.5). We 
then get the following exact formula for 4jo{x): 

See, for example, Montgomery and Vaughan [21, p. 401]. Further, Davenport [7, Equation 11, 
p. 110] gives the following bound for the error as a function of x and T: 

\tk\<T ^ k=l 

where x is an integer and is the imaginary part of pk, the k*^^ complex zeta zero. Note that, 
for a fixed x, the error term approaches as T approaches oo, consistent with (8.6). 

The second sum in (8.6) is too small to have any visible effect on the graphs that we will 
display, so we will use M = when we graph this formula. 

The dominant term, x, shows roughly how fast ijj{x) increases. In fact, it is known that iplx) 
is asymptotic to x (written ^p{x) ~ x), which means that lim^_j.oo = 1- 

The first two residues, x and arising from poles at s = and s = 1, give a good linear 

approximation to ^(x). That approximation, which comes from equation (8.6) with = and 
M = 0, is shown at the top left of Figure 6. This smooth approximation, while good on a large 
scale, misses the details, namely, the jumps at the primes and the prime powers. To replicate 
those details, the approximation needs some of the terms with zeta zeros in equation (8.6). 

Ingham [14, pp. 77-80] also proves equation (8.6). On page 80, he makes the intriguing 
remark. 

The 'explicit formula' [(8.6)] suggests that there are connections between the 
numbers (the discontinuities of ipolx)) and the numbers p. But no relationship 
essentially more explicit than [(8.6)] has ever been established between these two 
sets of numbers. 

This is still true today. 

It is also worth pointing out that in (8.6), the sum over complex zeta zeros does not require 
the calculation of the zeta function itself. In the approximations in later sections of this paper. 
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the sums do involve zeta or its derivatives, which makes it harder to prove the existence of an 
exact formula like (8.6). 

8.1. Numerical Experiments. We will now work through some numerical examples. Given 
a rectangle, we will compute line integrals on each of its sides. We will see that the integral on 
the first segment, from c — iT to c + iT, is significant, and that the integrals on the other three 
sides are small. This means that the first integral alone is approximately equal to the integral 
around the entire rectangle, which, in turn, equals the sum of the residues inside the rectangle. 
(Note that Perron's formula, equation (8.2), includes a factor of l/(27ri), so the right side of 
(8.2) is the sum of the residues, not 2m times the sum of the residues.) As we integrate around 
the rectangle in a counterclockwise direction, we'll let Ji, I2, I3, and I4 be the line integrals on 
the right, top, left, and bottom sides, respectively. 

The Dirichlet series (8.1) converges for 3fJ(s) > 1, so in the integral in (8.2), we must have 
c > 1. Let's take c = 2. Let's extend the left-hand side of the rectangle to 3ft(s) = —1, the 
rectangle encloses no real zeros. Then we will have M = in (8.6). 

Equation (8.2) suggests that, for a given x, if T is large, then the sum of the residues may be 
close to ipo{x). We will experiment with different values of T in (8.2), that is, different values 
of in equation (8.6). 

Let's start with x = 10. Then ^jq^x) = ij{x) = log(2520) ^ 7.832. If T = 14, then the 
lower-right corner of the rectangle is at 2 — lAi, and its top left corner is at — 1 + lAi. The first 
complex zeta zero is at s ~ 1/2 + 14.135z, so this rectangle encloses no complex zeta zeros. 
Therefore, we will have = in equation (8.6). Numerical integration tells us that Ji ~ 8.304, 
I2 ^ -.069 - 0.305i, I3 ^ -0.004, and h ~ --069 + 0.305i. The sum of these four integrals is 

/l+/2 + /3 + /4^ 8.162. 

As a check on our numerical integrations, we can substitute x = 10, A^ = 0, and M = into 
equation (8.6). This, in effect, computes the sum of the four integrals by using residues. We get 
8.162, as expected. Cauchy was right: we can compute a contour integral by adding residues! 

Now let's extend the rectangle in the vertical direction to enclose one pair of complex zeta 
zeros. We will now have A^ = 1 in equation (8.6). We can, for example, take T = 15, because the 
first two zeta zeros are 1/2 + 14.135i and 1/2-1- 20.022i. A computation of the integral over the 
sides of this rectangle gives h ^ 7.602, I2 ^ .078 - 0.213z, I3 ^ -0.006, and h ~ -078 + 0.213i. 
The total is ~ 7.751, a little closer to our goal, ipo{x) ~ 7.832. 

Because the second zeta zero has imaginary part roughly 20.022, we could have enclosed one 
pair of zeros by extending our rectangle up to T = 20 instead of to T = 15. Had we done this, 
we would have found that h ^ 7.518, h ^ .119 + 0.087z, I3 ^ -0.005, and h ~ -119 - 0.087i 
However, the sum of these four integrals is unchanged (~ 7.751), because the rectangle with 
T = 20 encloses exactly the same poles as the rectangle with T = 15. 

Let's set T = 100 and extend the rectangle again, so its lower-right and upper-left corners 
are 2 — lOOi and — 1 -|- lOOi. There are now A^ = 29 zeta zeros with positive imaginary part 
less than T. The integrals over the sides of this rectangle are Ii ~ 7.815, I2 ~ —.014 — 0.018i, 
J3 ^ -0.005, and h ~ --014 + 0.018^. The sum is ^ 7.782, not far from the goal of 7.832. 
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In Figure 6 below, we show the step function ip{x), along with approximations to tplx) using 
formula (8.6) with M = 0. To create the four graphs, we set = 0, = 1, = 29, and 

= 100. We computed (8.6) for x = 2 through x = 50. Notice that, if x is a value that causes 
the step function to jump, the smooth curve is usually quite close to the midpoint of the y 
values, that is, ipo{x). 



N=l 
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Figure 6. approximating ip{x) with A^ = 0, A^ = 1, A^ = 29, and A^ = 100 pairs 
of zeta zeros 

In Figure 6, notice that for a fixed A^ (or T), the agreement gets worse as x increases. Second, 
for a given x, as A^ (or T) increases, we will usually see better agreement between the plotted 
points and the step function. 

What would happen if we used this integral approximation, 

^o(x)^— / --^-ds, (8.8) 

based on (8.2), instead of the summing residues? In Figure 7, we set c = 2 and T = 100 in (8.8) 
and perform the integration for many x < 50. Recall that a rectangle with T = 100 encloses 
A^ = 29 pairs of complex zeta zeros, so in this sense, it is fair to compare this approximation 
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with the third graph in Figure 6. We can see that, at least in this example and for these x, the 
approximations are about equally good. 




Figure 7. Approximating iIj{x) with an integral: c = 2, T = 100 

Incidentally, there's nothing special about c = 2 here. If we had computed the integral using 
c = 1.1 or c = 1.5 instead, the graphs at this scale would look about the same. If we had used 
c = 2.5, the graph would be a little bumpier. 




Figure 8. x = 10.5, c = 1.1 

Here is another way to visualize how close the Perron integral (8.8) is to the summatory 
function. Take x = 10.5. ipoi^) — 7.832. Set c = 1.1 in equation (8.8). Figure 8 shows how 
close the integral is to 7.832 as a function of T. 

Again taking x = 10.5, figure 9 shows the integrals with c = 1.1 and c = 2.0. 

Notice that the graphs in Figure 6 do not show any sign of Gibbs phenomenon. However, 
this is due only to the scale of the graphs. Figure 10 shows graphs similar to those of Figure 6, 
but with = 100 and A^ = 400, and for x between 15 and 20. The Gibbs phenomenon at 
X = 17 is quite clear: the step function jumps by '?/'(17.1) — '?/'(16.9) = log(17) ~ 2.833. At 
X = 16, the Gibbs phenomenon is less discernible because the jump in the step function there 
is only V^(16.1) - V(15.9) = log(2) ^ .693. 
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Figure 9. x = 10.5 : c = 1.1 (blue) and c = 2.0 (red) 



As with Fourier series, the undershoot and overshoot seem to be related to the size of the 
jump. For example, let's see what happens at x = 16, x = 17, and x = 97 and x = 127, 
where the step function ip^x) has jumps of size log(2), log(17), log(97) and log(127). We'll use 
equation (8.6) with = 2000 pairs of zeta zeros. For each of these values of x, let i/min be 
the y-value at the minimum just before the jump, and let i/max be the y-value at the peak just 
after the jump. Taking into account the size of the jump, the number that corresponds to the 
Gibbs constant g in equation (5.1) is: 

7 Umax Umin C 1 



where c is the size of the jump. (The final division by c normalizes for the size of the jump.) 
For X = 16, x = 17, X = 97, and x = 127, d is about 0.08965, 0.08951, 0.08957, and 0.08974. 
2000 terms are not enough to see the small jump of size log(2) at x = 128. 




Figure 10. The Gibbs phenomenon in approximations to tplx) 
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9. Computing the Mertens Function M{x) 

Recall the definition (4.3) of the Mobius /i function. The summatory function of fi{n) is the 
Mertens function, usually denoted by M{x): 

X 

M(x) = 5^/i(n). 

n=l 

The Mertens function has an interesting connection to the Riemann Hypothesis. In 1897, 
Mertens conjectured that |M(x)| < -y/x for all x > 1. It is known [17] that the Mertens 
conjecture is true for x < 10^^. If this conjecture were true, then the Riemann Hypothesis 
would also have to be true. But in 1985, Odlyzko and te Riele proved that the Mertens 
conjecture is false [22]. (This does not disprove the Riemann Hypothesis.) Odlyzko and te 
Riele did not give a specific x for which the Mertens conjecture is false. In 2006, Kotnik and 
te Riele proved [17] that there is some x < 1.59 x lO^'' for which the Mertens conjecture is 
false. The 1985 results used sums with 2000 pairs of complex zeta zeros, each computed to 100 
decimals. The 2006 results used sums with 10000 pairs of zeros, with each zero computed to 
250 decimals. 

Here is the Dirichlet series involving /i(n): 



E 

n=l 



This holds for ^{s) > 1. 

From this Dirichlet series. Perron's formula suggests an approximation to the summatory 
function. As usual, the prime in the summation and the subscript below means that, where 
the step function M{x) has a jump discontinuity, Mq{x) is the average of the function values 
just before and after the jump. 

y u(n) = Mo(x) = lim / -——ds. (9.2) 

We will approximate this line integral with an integral over a rectangle that encloses poles of 
the integrand in (9.2). We will first locate the poles of the integrand, then compute the residues 
at those poles. The poles occur at the places where the denominator is 0, that is, s = 0, and 
at the zeros of zeta. The residue at s = is ^ = —2. 

To get the residue at pk, the k^^ complex zeta zero, we will use (6.1) for the residue of the 
quotient A{s)/B{s) where 

A{s) 1 x' 



B{s) as) s 

To use (6.1), we will set A{s) = ^ and B{s) = C{s). Then B'{s) = C{s), and the residue of 
ferf at s = Pk is 



((s) s 



X- 



Pk 



pkC {pk) ' 
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At the real zeros s = —2, —4, —6 . . . , we get a similar expression, but with replaced by 



The alleged approximation that uses the sum of the residues over the first complex pairs 
of zeros and the first M real zeros, is 



The second sum is small for large x. As we will see, only about 3 terms of this series are 
large enough to affect the graph, and then, only for small x. In Section 9.2 will discuss to what 
extent this approximation can be proved. 

Note that we have not proved that the expression on the right of (9.3) is close to Mq{x). 

We will form a rectangle that encloses the pole at s = 0, along with the poles at the first 
A^ complex zeros, and the poles at the first M real zeros. We will then integrate around the 
rectangle. We will see that the line integrals on the top, left, and bottom sides are small, so 
that the line integral in (9.2) can be approximated by the sum of integrals on all four sides, 
which, in turn, equals the sum of the residues of the poles that lie inside the rectangle. 

9.1. Numerical Experiments. In the first three graphs in Figure 11, we set M = 0, and 

used A^ = 0, A^ = 29, and A^ = 100 pairs of complex zeta zeros in (9.3). Notice that, even 
with A^ = 100, the approximation is not very good for small x. So, in the fourth graph, we 
also include the sum of the residues of the first M = 3 pairs of real zeta zeros. This makes the 
agreement quite good, even for small x. 

Now let's look at numerical examples. Let's set c = 2, and use (9.3) with A^ = 29 and M = 3. 
A rectangle that encloses the pole at s = 0, the poles at the three real zeros s = —2, s = —4, 
and s = —6, and the 29 pairs of complex zeros could extend, for example, from 2 — lOOz to 
—7 + lOOi. If X = 32, then one can see from the graph that M{x) = —4, and, since this x is 
divisible by a prime power, /i(x) = 0, so that M{x) is continuous at this x. The integrals along 
the right, top, left, and bottom sides are, respectively, Ji ~ —4.506, I2 ~ .180 — 0.325z, ~ 0, 
and h ~ .180 + 0.325z. The total is ^ -4.145. 

Notice that, for a fixed A^ and M, (9.3) diverges more and more from the step function as x 
increases. We can remedy this by using a larger value of A^ in (9.3). 

In Figure 11, we drew approximations based on the sum of the residues given by equation 
(9.3). What would happen if we chose c and T, then, for each x, we plotted the integral 



instead of the sum of the residues? 

On the left side of Figure 12, we show the values we get by substituting c = 2 and T = 100 in 
(9.4). The integral (9.4) is clearly following the step function, but the agreement with the step 
function is not very good. On the right side of Figure 12, we plot the integral using c = 1.5. 
With c = 1.5, the agreement is much better. In fact, for the x in range of this graph, this 
integral approximation is about as good as the sum of residues when A^ = 29 pairs of zeros are 
included. 



-2k. 




(9.3) 




(9.4) 
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Figure 11. approximating Mq{x) with = 0, = 29, and = 100 pairs of 
zeta zeros 



c = 2, T = 100 



c= 1.5, T= 100 




Figure 12. approximating Mo{x) with integrals: T = 100, c = 2 and c = 1.5 

Since the integral approximation is so bad with c = 2, T = 100, one might ask whether we 
can get a better approximation by taking a larger T. So, let's try a T that corresponds to, say, 
100 pairs of zeros; to do this, we can set T = 237. The graph is in figure 13. The overall level 
of agreement isn't that different, but using a larger T introduces higher-frequency terms into 
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c = 2, T = 237 



c = 1.5, T = 237 




Figure 13. approximating Mq{x) with integrals: T = 237, c = 2 and c = 1.5 

the graphs. What does improve the approximation is to use smaller a value of c (1.5 instead of 
2). 

N = 100, M = 3 

01 — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ^ — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' 

13 14 15 16 17 




Figure 14. Gibbs phenomenon in M{x): N = 100, M = 3 

Figure 14 shows the Gibbs phenomenon in the graph of M{x). 

9.2. Theorems About Mq{x). It certainly appears from the above graphs that the sums in 
(9.3) provide a good approximation to the step function Mq{x). Here we summarize what has 
actually been proven. 

Assuming the Riemann Hypothesis, and assuming that all zeta zeros are simple (i.e., are of 
multiplicity 1), Titchmarsh [33, p. 318] proved the following Theorem: 



Theorem 1. There is a sequence T^, v <Tu < v + 1, such that 



Mq{x) = -2 + lim V ^— + V 



\f-'^(2Txlx) 



2k 



h\<T. 



pCip) ^ (2fc)!fcC(2A: + l) • 



(9.5) 
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Note: because 



C'(-2A:) = (-1)^-|^C(2A: + 1), 



2(27r)2fc 

the /c*'* term in the second sum is the same as the k*^^ term in equation (9.3), namely, 



[-2k)c'i-2ky 

In Titchmarsh's Theorem, 7 is the imaginary part of the zeta zero p. The proof in [33] sums 
the residues inside a rectangle, the sides of which must not pass through any zeta zeros. This 
means that the in this Theorem cannot be the imaginary part of a zeta zero. In addition, 
the proof makes use of a theorem which says that, given e > 0, each interval [u, u + 1] contains 
a T = Tj, such that 

^ O(T^) 



C(a + zT) 



for I < (T < 2. 



Bartz [1] proves a result similar to Titchmarsh's Theorem, but without assuming the Riemann 
Hypothesis. Her result is also more general, in that it works even if the zeta zeros are not simple. 
In the easiest case, where the zeros are simple, she proves that there is an absolute positive 
constant Tq and a sequence T^,, T'^^Tq <Ty < (for u > 1), such that (9.5) holds. 

One would like to have this result 



1 ^ 1 



,=1 PkC'iPk) ^ i-2k)C'i-2ky 
but note that neither Titchmarsh's nor Bartz's result is this strong. 

10. Counting the Squarefree Integers 

A squarefree number is a positive integer that is not divisible by the square of a prime. The 
first 10 squarefree numbers are 1, 2, 3, 5, 6, 7, 10, 11, 13, and 14. 

The number of squarefree integers from 1 through x is usually denoted by Q{x). For example, 
Q(10) = 7, because 4, 8, and 9 are divisible by squares of primes, so they don't contributes to 
the count. The graph of Q{x) is an irregular step function. 

Define q{n) = 1 if n is squarefree and otherwise {q{l) = 1). Then Q{x) is the summatory 
function of q{n). 

The following Dirichlet series holds for 3f?(s) > 1 [13, p. 255]: 



E 

Perron's formula tells us that, if c > 1, 



g(n) C{s) 



..1 C(23) 



E' 



n=l 



q{n) = Q,{x) = hm — / ds (10.1) 

T^oo 2711 J^_iT C{2S) S 
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where Qo{x) equals Q{x) except where Q{x) has a jump discontinuity (that is, at the squarefree 
integers), and, at those jumps, Qo{x) is the average of the values before and after the jump. 

The integrand has a pole at s = 0, where the residue is 1. Recall that zeta has a pole of 
order 1 at s = 1. Therefore, the integrand has a pole at s = 1. The residue there is Gx/vr^. 

Finally, the integrand has poles wherever 2s is a zero of zeta. If p is any zeta zero, then the 
pole occurs at s = p/2. To get the residue at s = p/2, we will use equation (6.1) for the residue 
of the quotient A{s)/B{s) 

A{s) ^ Cis) X' 
B{s) C(2s) s ' 

where A{s) = C{s)^ and B{s) = C{2s). Then B'{s) = 2({2s), and the residue at the zeta zero 
s = p/2 is 

A{s) C(f) C(f) 



B'is) 2C'(p) f p C'ipY 
This holds whether p is a real or a complex zeta zero. Each zeta zero gives rise to one term 
of this form. Therefore, the sum over the residues at the first N complex zeros and the first M 
real zeros, plus the residues at s = and s = 1, gives this alleged approximation 

In the second sum, only the terms with M < 3 are large enough to affect the graph. Moreover, 
if k is even, then the term in the second sum 



_ 2fc 
X 2 



(^) x-'^a-k) 



{-2k)C'i-2k) {-2k)Ci-2k) 

is zero, because C(~^) = since —k is a zeta zero if k is even. Therefore, the second sum 
reduces to a sum over odd k < M. 

If the zeta zeros are simple, Bartz [2] proves that there is an absolute positive constant Tq 
and a sequence T^, 2^~^Tq <Ty< 2^To, (for z/ > 1), such that 

Qo(x) = 1 + ^ + hm y x^^i|^ + ya;-f-4i^p— , (10.3) 
^' ^^^^l^l^. P'^^P'^ & (-2fc)C'(-2fc)' ^ ' 

where 7 is the imaginary part of the zeta zero p. Bartz has a similar result [3] for cube-free 
integers. 

10.1. Numerical Experiments. Figure 15 shows graphs of two approximations to Qo{x) 
using (10.2). The graph on the left shows the linear approximation that comes from the first 
two terms of (10.2). The graph on the right shows (10.2) with N = 29, M = 3. The three 
terms in the second sum have a noticeable affect on the graph only for x < 10. 
The left side of Figure 16 shows a graph of the integral approximation 



2711 J^_iT C(2s) s 
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Figure 15. Approximating Qo{x) with sums of residues 
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c= 1.5, T= 100 
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Figure 16. Approximating Qo{x) with integrals; c = 2 and c = 1.5 

with c = 2 and T = 100. As noted above, a rectangle with T = 100 will enclose 29 pairs of 
complex zeros. The right side shows the integral approximation, this time, with c = 1.5. 



11. Tallying the Euler Function 

The Euler phi (or "totient") function is defined as follows. If is a positive integer, then 
(f){n) is the number of integers from 1 through n that are relatively prime to n. In other words, 
0(n) is the number of integers k from 1 through n such that the greatest common divisor of k 
and n is 1. For example, 0(10) = 4, and if n is prime, then 0(n) = n — 1. We will let ^(x) 
denote the summatory function of 0(n), so that 



X 

n=l 
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Then 

X 

n=l 

SO that $0(2^) equals except where has a jump discontinuity (this happens at every 
positive integer since 0(n) > 0). Further, at those jumps, ^o{x) is the average of the values of 
$(x) before and after the jump. 

The following Dirichlet series holds for 3?(s) > 2 [13, Theorem 287, p. 250]: 

n=l ^ ' 

From this Dirichlet series. Perron's formula tells us that, if c > 2, then 

The integrand has a pole at s = where the residue is 1/6. C,{^s) has a pole of order 1 at 
s = 1, so, because of the C(s — 1) in the numerator, the integrand has a pole at s — 1 = 1, that 
is, at s = 2. The residue at s = 2 is Sx^/tt^. 

The integrand also has poles at every real and complex zero of zeta. We will apply residue 
formula (6.1) with B{s) = ({s) and A{s) equal to everything else in the integrand, that is, 
A{s) = ({s — l)x^/s. The residue at the zeta zero s = p will be 

A{p) _C{p-l)xP 1 _x%{p-l) 
~ P W) ~ pC'(p) 

Therefore, if we create a rectangle that encloses the first N pairs of complex zeta zeros, the 
first M real zeros, and the poles at s = and s = 2, then the integral around the rectangle is 
just the sum of the residues, which suggests that 

A quick calculation will show that, for x > 1.5, the second sum is too small to visibly affect 
our graphs, so we will take M = in the graphs below. 

Notice that, if no zeta zeros are included in (11.3) (that is, if = M = 0), then the estimate 
for $(x) is given by the quadratic 

6 TT^ 

Figure 17 shows this quadratic approximation to $(x). This approximation looks quite good. 
Not only that, in Figure 17, the value of this quadratic appears to be between $(a; — 1) and 
$(a;) for every integer x. In fact, the first integer x where 

Hx-l)<l + ^<^x) 
6 7r^ 
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Figure 17. Approximating with the quadratic | + ^ 



is false, is x = 820: $(819) = 204056, $(820) = 204376, but | + ~ 204385.25831. There 
is only one such x less than 10'^, but there are 36 of them less than 10^, 354 of them less than 
10^, 3733 of them less than 10^, and 36610 such x less than 10^. 

Compare this to the standard estimate with an error term in [13, Theorem 330, p. 268] that 



E3x^ 
(pin) = — — + 0(x\ogx). 



n=l 



If (11.3) is a valid approximation to ^o{x) (or $(x)), (11-4) implies that 



\ PkC {pk) 



0{x log a;) 



11.4) 
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N = 29, M = 
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c = 2.1,T= 100 
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Figure 18. Approximating ^o{x) with a sum and an integral 
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Here is the integral approximation that we get from equation (11.2): 

The graphs in Figure 18 show that the sum (11.3) and the integral (11.5) seem to give 
reasonable approximations to the step function ^{x) (or ^o{x)), at least for these x. Quite a 
bit of analysis would be required to estimate how close (11.3) and (11.5) are to ^o{x)- 

R§kos [26] discusses this explicit formula related to Euler's (p function: 



f{z) = lim V 



n-!>oo ^ ' C'(p) 

0<S{p)<T„ ^ 

12. Tallying the Liouville A Function 

The Liouville lambda function is defined as follows. For a positive integer n, factor n into 
primes. If fl{n) is the number of such factors counting multiplicity, then A(n) = (— 1)*^^"^ For 
example, if n = 12, then X{n) = (—1)^ = —1, because 12 has three prime factors: 2, 2, and 3. 
We will denote the summatory function of A by L{x): 

X 

L{x) = J2Hn). 

n=l 

It is not hard to see that L{x) is the number of integers in the range 1 < n < x that have 
an even number of prime factors minus the number of integers in that range that have an 
odd number of prime factors. This is because each n with an even number of prime factors 
contributes +1 to the sum L{x), while each n with an odd number of prime factors contribute 
-1. 

The following Dirichlet series holds for 3fJ(s) > 1 [13, Theorem 300, p. 255], [9, eq. D-53], 
[20, p. 255]: 

J2 ^(^) - C(2^) 



The analysis is similar to those in the previous sections. 

From this Dirichlet series. Perron's formula tells us that, if c > 1, then 

Lo(x)= lim — r^^^^—ds. (12.1) 

The integrand has a pole at s = where the residue is 1. C,{s) has a pole of order 1 at 
s = 1, so, because of the C(2s) in the numerator, the integrand has a pole at 2s = 1, that is, at 
s = 1/2. The residue at s = 1/2 is y^/C (l) ^ -0.684765v^. 

The integrand also has poles at every real and complex zero of zeta. We will apply residue 
formula (6.1) with B{s) = ({s) and A{s) equal to everything else in the integrand, that is, 
A{s) = ({2s)x'^/s. The residue at the zeta zero s = p will be 

A{p) ^ C{2p)xP 1 ^ x^C(2p) 

B'{p) p C'(p) pC'(p) ■ 
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This expression holds for any zeta zero, whether real or complex. Recall that the negative 
even integers are the real zeros of zeta. Therefore, if p = —2k is any of these real zeros of zeta, 
then the C(2p) in the numerator is always 0, because 2p is another real zero. So, in the sum of 
the residues that correspond to the real zeros, every term will be 0. 

Therefore, we can create a rectangle that encloses the first N pairs of complex zeta zeros and 
the poles at s = and s = 1/2. When we integrate around the rectangle, the integral is just 
the sum of the residues, which suggests that we might have the approximation 



Lo(x) 



X 



C(i) 



N 



.fc=i 



. C(2pfc) 

pkC (pk) 



(12.2) 



N = 



N = 




Figure 19. Approximating Lo{x) with 1 + -^^^ 

Figure 19 shows the approximation to the step function based on the first two terms of 
equation (12.2). As you can see, (unlike the corresponding approximation in the previous 
section), these first two terms completely obscure the behavior of the step function. 



N = 29 



N= 100 




Figure 20. Approximating Lq{x) with sums over iV = 29 and N = 100 pairs of 
zeta zeros 



30 



ROBERT BAILLIE 
c = 1.5, T= 100 




Figure 21. Approximating Lo{x) with an integral; c = 1.5, T = 100 



N = 29 




200 



-10 - 



-15 - 



Figure 22. Approximating Lq{x) with a sum over = 29 pairs of zeta zeros 



You can see from the graph on the left of Figure 19 that for 1 < x < 50, the step function L{x) 
is less than or equal to 0. In the graph on the right, L(96) = L(586) = 0, so for 1 < x < 1000, 
L(x) is also less than or equal to 0. In 1919, George Polya conjectured that L{x) < for all 
X > 1. However, in 1958, this conjecture was proven to be false. We now know [31] that the 
smallest counterexample x > 1 for which L{x) > is x = 906150257. The first ten values of x 
for which L(x) = are 2, 4, 6, 10, 16, 26, 40, 96, 586, and 906150256 [30, Seq. A028488]. It is 
not known [31] whether L{x) changes sign infinitely often, but it is known (see [31] and [32]) 
that 



[X 



X 



does change sign infinitely often. 

For small x, the approximations in Figure 20, are not terribly close to the step function 
Lo(x), even when 100 pairs of zeros are used. However, the integral approximation shown in 
Figure 21 is a lot better. 
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To investigate this discrepancy, let's set x = 3.5. The value of the step function at a; = 3.5 is 
L(3.5) = —1. Let's integrate around the rectangle that extends from 1.5 — lOOz to — 1 + lOOi. The 
integrals along the right, top, left, and bottom sides are: Ji ~ —1.01558, 12 ~ 0.00623— 0.01804i, 
J3 ^ .35419, and h ~ 0.00623 + 0.01804i. The sum of these four integrals is ^ -0.64893. Note 
that Ji, the integral from 1.5 — lOOi to 1.5 + lOOi, is fairly close to L(3.5). If the integral along 
the left side (J3) was, say, 0, then the sum would be ~ —1.00312, much closer to L(3.5) 

Finally, Figure 22 shows that, for x < 200, the sum over = 29 pairs of zeros tracks the 
irregular step function L{x) pretty closely. As usual, as x increases, for a fixed A^, the fine 
details of the step function gradually get lost. 



13. Tallying the Squarefree Divisors 

For positive integers n, v{n) is defined to be the number of distinct prime factors of n. 
Therefore, if the prime factorization of n is 

then u{n) = k. It is not hard to see that 2*^*^"^ is the number of squarefree divisors of n. This 
example shows why: Let n = 60 = 2^ ■ 3 ■ 5. i>{60) = 3, because 60 has three distinct prime 
factors. We can list all the squarefree divisors of 60 by taking all 2^ = 8 combinations of the 
presence or absence of the three primes, that is, 2" ■ 3^ ■ 5^^, where a, b, and c independently take 
the values and 1. 

The following Dirichlet series holds for > 1 [13, Theorem 301, p. 255], [9, eq. D-9]. 



..1 C(23)- 
Let's denote the summatory function of 2*^*^") by T(x), so that 



n=l 

Perron's formula says that, for c > 1, 



1 r^'+^'r c(s)'^ x^ 
To(x)= lim — / ^^—ds. (13.1) 

where the subscript has the usual meaning. The integrand 

gsf 

C(2s) s 

has poles at s = and at every place where 2s is a zeta zero. It also has a pole of order 2 at 
s = 1, due to the presence of zeta squared in the numerator. The residue at s = is —1/2. 
The residue at s = 1 is calculated as follows. 

According to (6.4), the residue can be obtained by differentiating the product of (s — 1)^ 
and the above integrand, then taking the limit as s approaches 1. When we differentiate the 
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product 



(s-1 



C(2s) s ' 
with respect to s, we get 

_{s- 1) VC(s)2 2{s - 1) VC(s)2C'(2s) {s - l)'^xX{sf log(x) 
s^C{2s) sC(2s)2 ^ sC(2s) 

2(s- l)VC(s)C'(s) 2(s- l)x^C(s)^ 
^ <(2s) ^ sC(2s) ■ 

The hmits of the first three terms as s approaches 1 are easy because [13, Theorem 281, p. 
247] 

hm(s - l)C(s) = 1. 

({2) = 7r^/6, so the hmits of the first three terms are — Gx/vr^, —72x(' {2) / tt'^ , and 6xlog(x)/7r^. 
For the fourth and fifth terms, we will need two facts. First [8, p. 166], 

(is) = , ^ + 7 + E S^7n(^ - 1)" (13.2) 

s — 1 ^-^ nl 

^ ' n=l 

where 7 ^ .57721566 is Euler's constant, and 7„ are the Stieltjes constants. Second [13, 
Theorem 283, p. 247], 

C'W = ^ + o(i), 

The sum of the 4th and 5th terms above is 

<(2s) 

The limit of 

2x-(g-l)C(g) 
<(2s) 

is 12x/7r^, and we claim that the limit of 

C(s) + (s-i)C'(.) 

is 7. To see this, we write out the first few terms of the Laurent series for C, and C,': 

as) + - i)c'(^) 

= + 7 + (powers of (s - 1)) + (s - 1) ( + 0(1' 

s — 1 \s — 1 

= 7 + (powers of (s - 1)) + (s - 1)0(1). 

As s approaches 1, this expression approaches 7, as claimed. Combining these limits, the 
residue of the integrand at s = 1 is 

-Qx 72xC'(2) 6a;log(a;) 12a;7 62; (-12C'(2) + 277r2 - tt^) 6a;log(a;) 
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Next, if p is any real or complex zeta zero, we can derive an expression for the residue at 
s = p/2 using the same method as in section 10 where we counted the squarefree integers. 
Again, we will use equation (6.1) for the residue of the quotient A{s)/B{s) 



B{s) C{2s) s ' 

where A{s) = C('S)^ v ^i^) — C(2s). Then B'{s) = 2({2s), and the residue at the zeta zero 
s = p/2 is 

A{s) C(f)'xf xHiir 



B'{s) 2C'{p) 



pC'(p) 



Putting all these residues together, we get this alleged approximation for T{x) that involves 
the first pairs of complex zeta zeros and the first M real zeros: 



N 



To{x) ~ ai + a2X + a3a;log(x) + 23? ( — 

* PkQ [Pk 



M 



_ 2fe 
X 2 . 



k=l 



-2k)C'i-2k) 



(13.3) 



where 



02 



ai = -1/2, 

6(-12C(2) + 277r^-7r^) 
^4 



0.78687, 



and 



as = 6/7r^ 



.60793. 



The second sum is too small to visibly affect the graphs, so we'll set M = when we draw 
the graphs below. 



N = 0, M = 




Figure 23. Approximating the sum of 2*^*^"^ with the first three terms of (13.3) 
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N = 100,M = c=1.5,T = 237 




45 50 55 60 45 50 55 60 

Figure 25. Approximating the sum of 2^^^") with the same sum and integral for 



40 < X < 60 

We will now use (13.3) to draw a graph of Tq{x). Figure 23 shows the approximation to the 
step function based on the first three terms of equation (13.3). Figure 24 shows the approxi- 
mations using (13.3) with = 100, M = 0, and the integral 

with c = 1.5 and T = 237. This value of T is large enough for the corresponding rectangle 
to enclose 100 pairs of complex zeros. The approximations are quite good in this range of x. 
Figure 25 shows the same approximations, but for the range 40 < x < 60. 

Mathematica can compute the residue at the second-order pole at s = 1, as follows: 

res = Residue [Zeta[s] ~2/Zeta [2 s] * x"s/s, {s, l}] 
Collect [res, Log[x] , Simplify]. 
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The second command separates out the term that contains log(a;). The result is 

6xlog(x) ^ 6x ((27 - l)7r2 - 12C'(2)) 
^2 ^4 

Assuming the zeta zeros are simple, Wiertelak [37, Theorem 4] proves that, for x > 1/4, 

for a certain sequence T„, where Oi, 02, and 03 are the values given above. Note: because 

C'(-2^) = (-l)'^C(2fc + l), 

the second sum in (13.4) has the same as form the second sum in equation (13.3). 

If we have the Dirichlet series for an arithmetical function, the procedure for deriving an 
approximation should be clear by now. So, from now on, we will give fewer details about the 
derivations of the residues except where new issues arise, such as poles of order 2, 3, or 4. In 
each section below, we'll use the same symbol, T{x), for the respective summatory function of 
that section, and Tq{x) for the summatory function modified in the usual way. 

14. Computing sigma Sums a{n'^) and a{ny 

(j{n) is the sum of the (positive) divisors of the positive integer n. Here we consider two 
sums involving the a function that can be approximated using Perron's formula. 

14.1. The Sum of cr(n^). This Dirichlet series holds for s > 3 [9, set /c = 1 in eq. D-51], [20, 
eq. 5.39, p. 237]: 

a{n^) C(5)C(s-l)C(^-2) 



n=l 



C(2s - 2) 



Note that there is a typo in [9, eq. D-51]: the denominator on the right side should be 
C(2(s — k)). Let Tq{x) denote the (modified) summatory function of cr(n^). Then, based on 
Perron's formula, we have this integral approximation 

2vrz Jc-iT C(2s - 2) s 

We will not derive an estimate that shows how close this integral is to the sum for given 
values of x and T. However, we will see that, for modest values of x and T, the integral follows 
the step function rather closely. 

The integrand has poles ats = 0, s = l, s = 2, and s = 3. There are also poles at those s 
such that 2s — 2 = p, where p is any zeta zero, that is, at each s = p/2 + 1. The sum of the 
residues, which we hope will be close to Tq{x), is 

To(x) ^ ai + a,x + a,x' + a,x' + 2^ f V ^"^^^i ^ ^1 wf/t " 1 (1^.2) 

2 (f + 1) C (Pfc) J 
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where 

1 

ai = — ~ 0.0208333, 
48 

1 

02 = -— ~ -0.083333, 
az = -\ = -0.25, 

and 

04 = ^^1^ ~ 0.608969. 

The terms oi, a2X, Osx^, and 04^^ are the residues of the integrand at s = 0, s = 1, s = 2, 
and s = 3, respectively. One would expect to have another sum similar to the one in (14.2), 
but over the real zeros instead of the complex zeros. However, this sum is always 0: If > 1, 
then if p = —2k is the k^^ real zeta zero, then at least one of p/2 + 1, p/2, or p/2 — 1 will also 
be a real zeta zero, which causes the k^^ term to drop out. 

The factor of 2 in the denominator of the terms in (14.2) arises when we compute the residue 
of A{s) / B{s)^ where 

A{s) = as)as-ms-2)- 

s 

and 

5(s) = C(2s-2). 

When we use residue rule (6.1), we compute the derivative of B{s) with respect to s. This 
produces the factor of 2 that appears in (14.2). 

N = 



10000 - 



5000 - 




Figure 26. Approximating the sum of cr(r2^) with the first four terms of (14.2) 

Figure 26 shows the polynomial approximation to the summatory function Tq{x) that comes 
from the first four terms of (14.2). 

Figure 27 shows approximations to Tq{x) using a sum (14.2) and an integral (14.1). Note 
that the approximating sum increases roughly as x^. Moreover, the increases from one value 
of X to another can be both quite large, and quite irregular. From x = 23 to x = 24, the step 
function changes from 7449 to 9100, an increase of 1651. From x = 24 to x = 25, the step 
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Figure 27. Approximating the sum of cr(n^) using a sum and an integral 



function changes from 9100 to 9881, an increase of 781. In spite of these rapid and irregular 
steps, the integral approximation tracks the step function quite well. 

14.2. The Sum of cr(ri)^. This Dirichlet series holds for s > 3 [13, set a = 6 = 1 in Theorem 
305, p. 256]: 



n=l 



(nr as)as-irc{s-2) 



C(2s-2) 



This Dirichlet series is similar to the one for the sum of cr(n^), except here, there is a pole of 
order 2 at s = 2. The integral approximation to the summatory function To(x) is 

^^^^^ asKis - ims - 2) 



nix) 



2-ni 



ds. 



14.3) 



C(2s - 2) s 

From the residues at s = 0, s = 1, s = 2, s = 3, and at each s = p/2 + 1, we get an expression 
that may approximate To(a;): 



N 



Tq{x) ~ ai + a2X + a^x'^ + a^x'^ log(x) + a^x^ + 23ft j x ^ 



-iC(f + i)C(f)\(f -1) 



where 



2(f + 1) C'(p.) 



;i4.4) 



ai 



576 
1 



0-2 



24 



i -0.001736, 
0.041667, 



«3 = ^f^^-47 + l-21og(27r; 



-0.765567, 



1 

'4 



-0.25, 
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and 



as 



503) 
6 



1.00171. 



The terms Oi and are the residues of the integrand at s = and s = 1, respectively. The 
residue at s = 2 gives rise to OsX^ and 04^^ log(x). The term a^x^ is the residue at s = 3. 



N = o 



25000 - 



15000 - 



10000 




Figure 28. Approximating the sum of (t(?t,)^ with the first five terms of (14.4) 



N = 100 



c = 3.5, T = 237 
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Figure 29. Approximating the sum of cr(n)^ using a sum and an integral 

Figure 29 shows the approximation to the step function up to x = 30 based on (14.4) with 

= 100 pairs of zeta zeros, and the integral (14.3) with c = 3.5 and T = 237. 

Figure 30 shows the sum in (14.4), but going only up to a; = 20. Because the scale is different, 
we can see that the sum is rather "wavy" for small x. This raises the question: what would 
happen if we used more zeta zeros in the sum? 

Figure 31, shows the same sum approximation as in the previous two Figures, but here we 
use A^ = 200 pairs of zeta zeros. With more zeros included in the sum, the "wavyness" is even 
more pronounced. One wonders whether the sum in (14.4) converges to to the step function as 
A^ approaches infinity. 
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N = 100 



8000 - 
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2000 - 




Figure 30. Approximating the sum of aln)"^ up to x = 20 using = 100 pairs of zeros 



N = 200 



25000 - 



20000 - 



5000 
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Figure 31. Approximating the sum of a{n)'^ using A^ = 200 pairs of zeros 



100 - 



-200 - 




Figure 32. The sum of the first 89 (sohd) and 90 (dashed) zeta terms of (14.4) 



Figure 32 shows a closeup of the sum of the first 89 and 90 zeta terms of (14.4) (that is, 
ignoring the initial terms involving Oi through 05). Between a; = 10 and x = 11, the sum 
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becomes more wavy when the 90*'^ term is added to the sum. At x = 10.33, the sum of the 
first 89 terms (before taking the real part) is about 43.497 + 11.589z. The 90*'^ term is about 
101.756 + 9.082i. The real part of the 89*^^ term is quite a bit larger than the real part of any 
of the previous terms, or of the sum of the first 89 terms. This unusually large term is what 
causes the waves to appear around x = 11. As to why the distortion appears in the form of 
waves, we'll discuss that in Section 26. 



15. Computing tau Sums r(n^) and r(n)^ 

r(n) is the number of (positive) divisors of the positive integer n. Here we consider two sums 
involving the r function that can be approximated using Perron's formula. 

15.1. The Sum of r(ri^). It can be shown that r(n^) is the number of ordered pairs of positive 
integers (i, j) such that lcm{i,i) = n. 

This Dirichlet series holds for s > 1 [9, set = in eq. D-51], [20, eq. 5.39, p. 237], [15, 
Equation 1.105, p. 35]: 

n=l ^ ' 

The integral approximation to the summatory function Tq{x) is 

The integrand has a pole of order 3 at s = 1, along with poles of order 1 at s = and at each 
s = p/2, where p ranges over all zeta zeros. The residue at s = is 1/4. The residue at s = 1 
is 

6 (tt^ (-371 + 1 - 37 + 37^) + 144C(2)2 - 12n' (("(2) - C(2) + 37C(2))) 



X- 



+ xlog{xf — , 



3 

where 7 is Euler's constant and 71 ^ —.072816 is the first Stieltjes constant. The 71 comes 
from equation (13.2), the series for the zeta function whose coefficients involve the Stieltjes 
constants. 

The sum over all the residues gives an expression that may approximate the summatory 
function: 

ro(x) ~ ai + a2X + a3xlog(x) +a4xlog^(x) + 23f^ ^//^ ] ^ ' 



X PkC'{Pk)) ^ (-2fc)C'(-2fc)' 

(15.2) 



where 

= 1/4 = .25, 



a2 



and 



EXPERIMENTS WITH ZETA ZEROS AND PERRON'S FORMULA 

6 (tt^ (-371 + 37^ - 37 + 1) + 144C'(2)2 - Utt^ {("{2) + 3^C{2) - C(2))) 
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6 



0.12226, 



03 = : ^ 1.13778, 



04 



3 



0.30396. 



The second sum is too small to have any visible effect on the graphs, so when we graph 
equation (15.2), we will set M = 0. 



N = 0, M = 




Figure 33. Approximating the sum of r(n^) with the first four terms of (15.2) 



N = 100 



c= 1.5, T = 237 




Figure 34. Approximating the sum of r(n^) using a sum and an integral 

The derivation of the expression for the residue at a third-order pole can be quite tedious. 
Fortunately, Mathematica can easily compute this residue for us: 

res = Residue [Zeta[s] ~3/Zeta [2 s] * x"s/s, {s, 1}] 
resTerms = Collect [res, -[Log[x] , Log[x]"2}-, Simplify]. 
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Then, resTerms [ [3] ] , resTerms [ [2] ] , and resTerms [ [1] ] are the terms with x, a;log(x), 
and a;log(x)^, respectively. 

Figure 33 shows the approximation we get from setting = M = in equation (15.2). 
Figure 34 shows the approximations using 100 pairs of zeta zeros in (15.2), and the integral in 
(15.1). 

15.2. The Sum of r(ra)^. This Dirichlet series holds for s > 1 [13, set a = 6 = in Theorem 
305, p. 256] 



E 



1 - C(2s)' 

n=l ^ ' 

This Dirichlet series is similar to the previous one, except that here, the pole at s = 1 is of 
order 4. The integral approximation to the summatory function To{x) is 

2711 J^_iT C(2s) s 

The residue at s = is —1/8. The sum of the residues at s = and s = 1, and at each 
s = p/2, suggests that this approximation might be valid: 

To(x) ~ ai + a2X + a^x log(x) + a^x log^(x) + a^x log^(x) (15.4) 

/ ^ /■ (Pk\ 4 \ A/ ^ / 2fc\4 

\h p^^'^p^)) h (-2^K'(-2fc) 

where 

ai = -1/8 = -.125, 

a2 = ^( - 47r^((3 - 127i)C'(2) + 2C(^)(2) - 3C"(2) + 187^C'(2) - 127(C'(2) - C"(2)) 
+ TT^ (^471 + 7(4 - 127i) + 272 + 47^ - 67^ - l) 

- 1728C'(2)3 + 1447r2c'(2)(2C"(2) + 47C'(2) - C'(2))) ~ 0.46032 

where 72 — —0.00969 is the second Stieltjes constant (see equation (13.2)), 

6 (tt^ (-471 + 67^ - 47 + 1) + 144C(2)^ - 127r^ (^^2) + 47C'(2) - C'(2))) 

^6 

3(47 - l)7r2 - 36C(2) 

a4 = 

and 



as = ^ — ^ ■ ■ ^^^^ ^ 0.82327, 

a4 = ~ 0.74434, 



1 

ag = — ~ 0.10132. 



7r2 



The residue at s = 1 accounts for the expressions involving 02, as, a4, and a^. The very 
complicated residue at s = 1 can be computed with the Mathematica commands 
res = Residue [Zeta [s] ~4/Zeta [2 s] * x"s/s, {s, 1}] 
resTerms = Collect [res, -[Log[x] , Log[x]"2, Log[x]~3}, Simplify]. 




Figure 35. Approximating the sum of r(n)^ with the first five terms of (15.4) 




Figure 36. Approximating the sum of r(n)^ using a sum and an integral 



Figure 36 shows approximations based on the sum (15.4) and the integral (15.3) approxima- 
tions. Figure 37 shows the sum approximation with = 100 and A^ = 200. Notice that this 
latter approximation is quite wavy, much like the one for (j(n)^ in figure 31. Figure 37 makes 
one wonder whether, for any given x, the sum even converges as A^ approaches infinity. 

In 1916, Ramanujan [25] stated that 

T{n)'^ = + a^x log(x) + a^^x log^(a;) + a^x log^(x) + 0{x^^'') 

n<x 

and gave the same values of and stated above; presumably, he also knew the values of 02 
and 03. A proof of this formula appeared in [38]; see also [24]. The error term can be reduced 
to 0(x2+^); see Ivic, [15, Equation 14.30, p. 394]. Therefore, if the right side of (15.4) does, 
indeed, represent the summatory function, then this error bound is, in effect, a bound on a sum 
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N = 100, M = 



N = 200, M = 



800 r 



600 



400 



200 




800 r 



600 



400 



200 




Figure 37. Approximating the sum of r(n)^ using = 100 and = 200 pairs 
of zeros 



over zeta zeros. That is, we would have 

/ oo 



a; 2 



PkC (pk 



0{x 



16. Tallying a{n)T{n) 
This Dirichlet series for a{n)T{n) holds for s > 2 [9, set a = 1, 6 = in eq. D-58]: 

■ ff(n)r(n) C(s)^C(s - 1)^ 



E 



n=l 



n° 



The integral approximation to the summatory function Tq{x) is 



2m 



c+iT 



ds. 



,.-,r C(2.-l) 

The integrand has a pole of order 1 at s = 0, and poles of order 2 at s = 1 and at s = 2. 
There is also a pole at each value of s where p = 2s — 1, where p runs over all zeta zeros. 
Therefore, these poles occur at each s = (p + l)/2. 

The sum of the residues, which we hope will be close to To(x), is 

Tq{x) ~ ai + a2X + a^x"^ + a^x^ log(a;) (16.2) 



M 



X 2 



2(^)C'(P.) 



-2k+l 



k=l 



)C'{-2k) 



The apparently spurious 2 in the denominators appears for the same reason as the 2 in (14.2). 
The second sum is too small to affect our graphs, so we will set M = when we graph (16.2). 
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Here, 



ai = -— ~ 0.020833, 
48 



02 = 2 = 0-5, 



as 



TT^ (24C(3)C(2) + (47C(3) - C(3) - 4^(3))) 
144C(3)2 



-0.17540, 



and 



TT 



a4 



1.12549. 



72C(3) 

ai and are the residues at s = and s = 1, respectively, a^x"^ + a^x"^ log(x) is the residue 



at s = 2. 



N = 0, M = 




Figure 38. Approximating the sum of a(n)r(n) with the first four terms of (16.2) 



N = 50 




c = 2.5, T = 144 



10 



20 



30 



40 



Figure 39. Approximating the sum of a{n)T{n) using a sum and an integral 



46 



ROBERT BAILLIE 



Figure 38 shows the graph of the sum approximation (16.2) without using any zeta zeros (i. 
e., with N = M = 0). Figure 39 shows the graphs of the sum approximation with 50 pairs of 
complex zeta zeros, and the integral approximation (16.1). 

17. Tallying \{n)r{n) 
This Dirichlet series for X{n)T{n) holds for s > 1 [9, set = in eq. D-47], [20, p. 225]: 

^ X{n)a{n) _ C{2sf 

71=1 3 \ / 

The integral approximation to the summatory function Tq{x) is 

This integrand has a pole of order 1 at s = 0, where the residue is 1. There is a pole of order 
2 where 2s = 1, that is, at s = 1/2. There is also a pole of order 2 at each s = p, where p runs 
over all zeta zeros. Residues at poles of order 2 in the numerator were discussed in Section 13. 
There, we computed the residue by taking a limit, using equation (6.4). Here, the residue at 
s = 1/2 can be derived in the same way. It is 

(-C(^)+27C(^)-C'(^)) , v^log(x) 



Here, for the first time, we encounter poles of order 2 in the denominator of the integrand. 
To compute the residues at s = p, we will use equation (6.2). Set B{s) = ({s)^, and set A{s) 
equal to everything else in the integrand, 

A{s) = a2sr- 

s 

so that the integrand is equal to A{s)/B{s). Then we substitute this A{s) and B{s) into 
equation (6.2). 

The result is a very complicated expression which we shall not write in full here. The 
denominator of this expression turns out to be 

3s'{as)C"{s) + C{sfY. (17.2) 

Remember that we want an expression for the residue at s = p, a zero of zeta. But notice that 
when we substitute s = p, every occurrence of ({s) in (17.2) becomes 0. The denominator then 
simplifies to 

SpXip)'. (17.3) 

We process the numerator in the same way. The numerator that results from substituting 
A{s) and B{s) into (6.2) is very complicated, but can be separated into two terms. The first 
term has a factor of log(x): 

3sxnogix)Ci2sY {Cis)C"{s) + C{sf). 
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Again, because ({s) = when we substitute s = p, this simphfies to 

3px''\og{x)C{2pn'{p)'. (17.4) 
The second term in the numerator has a factor of but without the log(a;): 
-xX{2s)i3C{sY{C{2s)~AsC{2s)) + 3sa2s)asK"{s) 

+ as)isa2sK^'\s) + 3C"(s)(C(2.) - 4<'(2.)))). 
Because ({s) will be when we substitute s = p, this simplifies to 

- x%i2p) {SCipr (C(2p) - 4<'(2p)) + 3pC(2p)C'(p)C"(p)) • (17.5) 

For the residue at any given zeta zero p, here's what we end up with. We get a quotient, 
the denominator of which is (17.3). The numerator has two terms, one of which, (17.4), has a 
factor of log(a;), the other of which, (17.5), does not. We can write the quotient for this zeta 
zero as 

Fi(p)x' + F,(p)xnog(x) 

w> ■ <^^-«' 

There is one of these expressions for each zeta zero. We sum (17.6) over all zeta zeros. When 
we combine this sum with the residues at s = and s = 1/2, we get the following expression, 
which may approximate the summatory function, To{x): 

Ux) ^ a, + a,^x + a3V^log(x) + 25J (^g j , (17.7) 

where 

ai = 1, 

a. -1.187104, 

and 

as = — ^ ~ 0.234452. 
2C(1)^ 

The functions Fi, F2, and G are 

F,is) = -C(2.) «(2s)C"(.) + C'(^) (C(25) - 4<(2.))) , 



F^is) = sa2syc{s 



and 



G{s) = sWsf. 

The following Mathematica code will perform the above calculations. This uses the residueFormula 
function from Section 6: 
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A[s] = Zeta[2 s] "2 * x"s/s; 
B[s] = Zeta[s] "2; 
expr = residueFormula[2, s] ; 
expr2 = Together [expr /. Zeta[s] 
numer = Collect [Numerator [expr2] 
numer [ [2] ] /x"s (* this 

numer [ [1] ] /(x"s Log[x]) 
Denominator [expr2] 



-> 0] ; 

Log[x] , Simplify] ; 

is Fl[s] *) 

(* this is F2[s] *) 



(* this is G[s] *) 

Note that (17.7) has no sum over the real zeta zeros. This is because, if s = p is one of the 
real zeros -2, -4, then the presence of 2s in both Fi and F2 guarantees that every term will 
be 0, because if s is one of these zeros, then so is 2s. 



N = o 



-10 : 
-15 - 




Figure 40. Approximating the sum of A(n)r(n) with the first three terms of (17.7) 



N = 50 



c = 1.5, T= 144 




Figure 41. Approximating the sum of X{n)T{n) using a sum and an integral 

Figure 40 shows the graphs of the sum approximation (17.7) without using any zeta zeros. 
Because 03 is positive, the expression consisting of the first three terms of (17.7) will eventually 
be positive. In fact, ai + a2y/x + a3y/x\og{x) > if x >~ 104.1. This suggest an overall 
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tendency for the summatory function T(x) to be positive, but otherwise, the first three terms 
bear no resemblance to the actual graph of T{x). Figure 41 shows the graphs of the sum 
approximation using 50 pairs of complex zeta zeros, along with the integral approximation 
(17.1). 

18. Tallying A(n)2'^(") 

Note: see Section 13 ("Tallying the Squarefree Divisors") for a discussion of the sum of 2'^^^\ 
This Dirichlet series for A(n)2'^(") holds for s > 1 [9, eq. D-21], [20, p. 227]: 

A(n)2^(") _ C(2g) 

n=l ' 

The integral approximation to the summatory function Tq{x) is 

To(x)^-^ P^^-rf.. (18.1) 
The sum over residues that may approximate the summatory function is 



where 
and 



02 = ~ 0.46890. 

car 



The functions Fi, F2, and G are: 

F,{s) = -sa2s)C{s) - Cism^s) - 2sC{2s)), 
F^is) = <(2.)C'(s), 

and 

G{s) = s^is)'. 

The following Mathematica code will perform these calculations: 

A[s] = Zeta[2 s] * x~s/s; 

B[s] = Zeta[s] "2; 

expr = residueFormula[2, s] ; 

expr2 = Together [expr /. Zeta[s] -> 0]; 

numer = Collect [Numerator [expr2] , Log[x] , Simplify]; 

numer[[2]]/x"s (* this is Fl [s] *) 

numer [[l]]/(x"s Log[x]) (* this is F2 [s] *) 

Denominator [expr2] (* this is G[s] *) 
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N = 50 



c = 1.5, T= 144 



15 - 



10 - 



-10 - 



15 - 




Figure 42. Approximating the sum of \(n)2''^^^ using a sum and an integral 

In Figure 42, notice that, for small x, the sum is not very close to T{x). 
Note: For this series, the corresponding sum over real zeros appears to diverge. Here is the 
sum of the first five terms: 

Ci-A) C(-8) CX-12) C(-16) C(-20) 

x2C'(-2)2 2xX'(-4)2 3x^('{-6y 4x^C'{-8y 5x^%'{-10y 
8.61152 65.2338 605.915 6409.28 73829.4 



x^ 



x^ 



X" 



x° 



X 



10 



In the previous examples, we have extended the rectangle around which we integrate far to 
the left, to enclose the real zeros of zeta. However, it is not necessary to do this. Aside from 
the discrepancy at small x, it seems sufficient to make the rectangle enclose only the complex 
zeros and the poles at s = and s = 1/2. 

19. Tallying \{n)T{n^) 
This Dirichlet series for X{n)T{n^) holds for s > 1 [9, eq. D-49], [20, p. 234]: 



n=l 



as] 



3 ■ 



The integral approximation to the summatory function To(x) is 



Tnix) 



c+iT 



C{2sfx' 



ds. 



2ni J^_iT C(s)^ s 
The sum over residues that may approximate the summatory function is 



(19.1) 



To(.) . a. + a.V-^ + asv/S log(.) + 2^ (f ^lM±lM]^^^^ , 

(19.2) 



where 
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02 = ^^^4 ~ 1.24413, 

and 

as = — ^ ~ -0.160544. 
The functions Fi, F2, F3, and G are: 

F,{s) = ?,s\{2sfC{sf + 2C'(3)^ (C(2s) (4s2c"(2.) + C(2s)) + AsW2sf - 4<(2s)C'(2s)) 
- <(2s)C'(s) {sa2s)C^'\s) - 3Cis) (C(2s) - 4<(2.))) , 



F,{s) = -<(2s)C'(s) {3sa2sK"is) + 2Cis) {C{2s) - 4<(2s))) , 



and 

G{s) = 2s^C{sf. 

The following Mathematica code will perform these calculations: 

A[s] = Zeta[2 s] "2 * x"s/s; 

B[s] = Zeta[s] "3; 

expr = residueFormula[3, s] ; 

expr2 = Together [expr /. Zeta[s] -> 0]; 

numer = Collect [Numerator [expr2] , {Log[x] , Log[x]"2}, Simplify]; 
numer[[3]]/x~s (* this is Fl [s] *) 

numer[[2]]/(x"s Log[x]) (* this is F2[s] *) 
numer [[l]]/(x"s Log[x]"2) (* this is F3[s] *) 
Denominator [expr2] (* this is G[s] *) 




Figure 43. Approximating the sum of A(n)r(?2^) using a sum and an integral 
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Because 03 in (19.2) is negative, the expression consisting of the first three terms of (19.2) 
will eventually be negative. In fact, ai + a2\fx + Osa/x log(a;) < if x >~ 1717.9. This suggest 
an overall tendency for the summatory function T(x) to be negative. Figure 43 shows the 
graphs of the sum approximation using 50 pairs of complex zeta zeros, along with the integral 
approximation (19.1). 

Note: For this series, the corresponding sum over real zeros appears to diverge. Here is the 
sum of the first five terms: 

2C(-4)^ C(-8)^ 2C(-12)^ C(-16)^ 2C(-20)^ 
a;2C'(-2)3 a;4C'(-4)3 3a;6C'(-6)3 1x^C^'{-%f 5a;i0('(-l0)3 

_ 4.51601 135.899 12996. 2.73295 x 10^ 1.03183 x 10^ 

nf>^ /y>4 /y»6 /y8 /ylO 

%AJ JU xAJ JU JU 

20. Tallying \{n)T{nf 
This Dirichlet series for \{n)T{nY holds for s > 1 [9, eq. D-50], [20, p. 234]: 

2 



A(n)r(n)2 C(2s 



00 

n=l ^ \ / 

The integral approximation to the summatory function Tq{x) is 

The sum over residues that may approximate the summatory function is 
To(x) ~ ai + a2\/x + a-i\/x\og{x) + ai\/x\og{xY 



I 2^ I ^^Pfe -^i(Pfe) + -^2(Pfc)log(a^) + i^3(Pfc)log(a^)^ + i^4(Pfc)log(a^) j ^202) 



,fe=i 



where 



ai = -2, 



+ ^ (0 (< (0 - (0) ^ (0'^ - (0^) - 

where 7 is Euler's constant, and 71 ~ —0.072816 is the first Stieltjes constant, 

-37C(i)+C(i)+2C'(^) 
ag = ^ -0.510158, 

and 

04 = — ^ ~ 0.027484. 
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The functions Fi, F2, F^, F4, and G are: 



F^{s) = -15s'a2sfC"isr - sC{s)\l2C"{s){a2s){4s\"{2s) + ((2.)) 

+ sC{2s)(sa2s)C^'\s) - 4C(^)(«)(C(2«) - 4<(2s)))) + 5s'a2s)C{s)C{s) 
(2<(25)C(^)(.) - 3C"mi2s) - 4<(25))) - 2C'{sri-12sC'i2s){2s\"i2s) + ((2^)) 
+ 12s'C'i2s)' + a2s){-8s'C^'\2s) + 12s'C"{2s) + 3C(2s))), 

F^is) = sC{s){15s\{2sYC"{sf + 6C'(s)^(C(2s)(4.\"(2s) + C(2.)) + As^i^sf 

- 4sC{2s)C{2s)) - 4sC{2s)C'{s){sC{2s)e\s) - 3C"(«)(C(2«) - 4<(2s)))), 

Fsis) = -3s'a2s)Cisf (2<(2s)C"(s) + C'(«) (C(2«) - 4<'(25))) , 



The following Mathematica code will perform these calculations: 

A[s] = Zeta[2 s] "2 * x'^s/s; 

B[s] = Zeta[s]~4; 

expr = residueFormula[4, s] ; 

expr2 = Together [expr /. Zeta[s] -> 0]; 

numer = Collect [Numerator [expr2] , {Log[x] , Log[x]''2, Log[x]''3}, Simplify]; 

numer[[4]]/x~s (* this is Fl [s] *) 

numer[[3]]/(x"s Log[x]) (* this is F2[s] *) 

numer[[2]]/(x''s Log[x]~2) (* this is F3[s] *) 

numer [[l]]/(x''s Log[x]'^3) (* this is F4[s] *) 

Denominator [expr2] (* this is G[s] *) 

Note: For this series, the corresponding sum over real zeros appears to diverge. Here is the 
sum of the first five terms: 



F4s) = s'a2sfCis)', 



and 



G{s) = 6sX{sy. 



4C(-4)-^ 



eX'(-4)^ 3x%'{-6)* 



X' 



C'(-16)3 4C'(-20)3 



2.36826 283.112 278745. 1.16534 x 10*^ 1.44207 x 10 
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21. Tallying X{n)a{n) 
This Dirichlet series for X{n)a{n) holds for s > 2 [9, set A; = 1 in eq. D-47], [20, p. 232]: 

A(n)a(n) _ C(2g)C(2g - 2) 
^ C(^)C(^^ - 1) 

The integral approximation to the summatory function Tq{x) is 

1 r^*^ C(2s)C(2s - 2) , 

To{x)^— ^-j-^, -^-ds. 21.1 

2vr« Jc-iT C(s)C(s - 1) s 

The integrand has poles at s = 1/2 and at s = 3/2. The pole at s = 1/2 gives rise to the 
term involving ai in equation (21.2) below. The pole at s = 3/2 gives rise to the term involving 

02. 

The integrand also has poles at each s = and each s = + 1. To compute the residues 
at these poles, we will proceed in two steps. First, for the residue at each s = pk, 'we apply 
equation (6.1) with 

, _ C(2g)C(2s - 2) X- 

^^'^ " C(^-i) T 

and 

B{s) = as). 

Equation equation (6.1) tells us that the residue is 

A{s) _ ^X(2s)C(25 - 2) 



B'{s) sCis - l)C'(s) 
When we evaluate this at s = p^, we get 

, C(2pfc)C(2pfc-2) 

""'Pkapk-iK'iPky 
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For the residue at each s = + 1, we take 

C(2s)C(2s-2) X' 



Ais) 



as) 

and 

B{s) = as-i). 

Equation equation (6.1) tells us that the residue is 

A{s) _ C(25)C(2s-2) 



B'is) sC{s)C'{s - 1 
When we evaluate this at s = + 1, we get 

„p.+i c(2(Pfe + i))c(2pfc: 



{Pk + mpk + ix'iPkY 

Therefore, the formula for the sum of the residues, which may approximate To{x), is 



Tnix) 

where 



(x) - aiv^ + a^x'/^ + 23? {x'"'F^ (p^) + x'-^^^F^ (pfc)) j , (21.2) 



ai = -— - — --- ~ -0.274495, 

12C(-^)Ca) 

a. - ^3^1(1) - -0.105029, 

^.)_C(2.)C(2.-2) 



and 



sCis - l)C{s) ' 
C(2s)C(25 + 2) 



(s + l)C(s + l)C'(s)- 

Note the presence of C{2s) in both Fi and F2. If s = — 2/c is a real zeta zero, then so is 
2s = —4k. Therefore, the sum over the real zeta zeros vanishes. 

The following Mathematica code will perform these calculations. Parts A and B calculate 
the residues at s = p^ and s = p^ + 1, respectively. 
(* part A. integrand has a pole at each zero of Zeta[s] *) 
(* set B[s] = Zeta[s] ; set A[s] = everything else in the integrand *) 
A[s] = Zeta[2 s] Zeta[2s - 2] / Zeta[s - 1] * x~s/s ; 
B[s] = Zeta[s] ; 

residueFormula[l , s] / x"s (* = Fl [s] *) 



(* part B. integrand has a pole at s = (1 + each zeta zero) *) 
(* set B[s] = Zeta[s - 1]; set A[s] = everything else *) 
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A[s] = Zeta[2 s] Zeta[2s - 2] / Zeta[s] * x"s/s ; 
B[s] = Zeta[s - 1] ; 

residueFormula[l , s]/x"s /. s -> s+1 (* = F2 [s] *) 



N = 50 



c = 2.5, T = 144 



200 - 



100 - 



-200 h 




Figure 45. Approximating the sum of A(ra)cr(ra) using a sum and an integral 

Figure 45 shows the summatory function of X{n)(y{n) and an approximation using the first 
50 pairs of zeta zeros. 




Figure 46. x = 10.5, c = 2.1 



Here is another way to visuahze how close the Perron integral (21.1) is to the summatory 
function. Take x = 10.5. The summatory function at a: = 10.5 is 15. Set c = 2.1 in equation 
(21.1). Figure 46 shows how close the integral is to 15 as a function of T. 

Again taking x = 10.5, figure 47 shows the integrals with c = 2.1 and c = 2.5. 



22. Tallying X{n)a{nf 
This Dirichlet series for A(n)cr(n)^ holds for s > 3 [9, set /i = /c = 1 in eq. D-46], [20, p. 232]: 
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Figure 47. x = 10.5 : c = 2.1 (blue) and c = 2.5 (red) 



>^ A(n)a(n)2 _ C(2s)C(2s - 2)C(2s - 4) 
^ ~ C(s)C(s-2)C(s-l)2 ' 



n=l 



The integral approximation to the summatory function Tq{x) is 
To(a:) ~ ^ C(2.)C(2.-2)C(2.-4) x- 

'°^^^-2vrJ_^ C(^)C(^-2)C(.-1)^ 3 ^ ^ 

The integrand has poles at s = 1/2, at s = 3/2, and at s = 5/2. These residues give rise to 
the terms involving ai, 02, and 03 in equation (22.2) below. 

The integrand also has first-order poles at each s = pk and at s = + 2. These residues at 
these poles are computed the same way we computed the residues in Section 21, above. These 
residues account for the terms 

x'"'Fi{pk)+xP'+^F2{pk) 

in equation (22.2). 

The integrand also has a second-order pole at each s = pk + I- To compute the residue at 
these poles, we use the standard formula (6.2) with 

C(2s)C(2s-2)C(2s-4) X' 



Ais] 



and 

This accounts for 
in the next equation. 



X 



as)as-2) 
B{s) = c{s-ir. 



GiPk) 
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Therefore, the sum over the residues, which we hope will approximate the summatory func- 
tion, is 

To{x) ~ aix^/^ + aaa;^/^ + aax^^^ 

+ 2^ (f^^'^^FM+x^^-^^FM+x'^^^^ '''^''^ + i^4(p.)log(x) \ ^22.2) 



where 



ai = ^ ~ -0.431757, 

1440C(-i)C(-D\a) 

a2 = ^^^^ ~ 0.02883, 

36C(-i)C(i)\(i) 



and 



^ C(3)C(5) ^ _0.018646. 

8C(i)C(i)\(i) 

Fi, F2, F3, F4, and G are given by 

_ C(25)C(25 - 4)C(25 - 2) 

^^'^ - 2)C(. - im«) ' 

C(2g)C(2g + 2)C(25 + 4) 
(s + 2)C(s + lK(s + 2)C'(s)' 



F3(s) = -(s + l)C(s - l)C(2s)C(s + l)C(2s - 2)C(2s + 2)C"(s) 

- C{s){{s + l)C(2s)C(s + l)C(2s - 2)C(2s + 2)C'(s - 1) 
+ C(s - l)((s + l)C(2s)C(2s - 2)C(2s + 2)C{s + 1) 

+ as + l)(C(2^s + 2)(C(2. - 2)(C(2.) - 2(. + 1)C'(2.)) 

- 2{s + l)C(2s)C'(2s - 2)) - 2{s + l)C(2s)C(2s - 2)C'(2(s + 1))))), 



F,is) = is + l)C(s - l)C(2s)C(^ + l)C(2s - 2)C(2s + 2)C'(s), 

and 

G{s) = is+iras-iras+ircisr. 

The following Mathematica code will perform these calculations. Parts A, B, and C calculate 
the residues at s = p^, s = + 2, and s — pk + 1, respectively. 
(* part A. integrand has a pole at each zero of Zeta[s] *) 
(* set B[s] = Zeta[s] ; set A[s] = everything else in the integrand *) 
A[s] = Zeta[2 s] Zeta[2 s - 2] Zeta[2 s - 4] / (Zeta[s - 2] Zeta[s - 1] "2) * x"s/s ; 
B[s] = Zeta[s] ; 

residueFormulaCl , s] / x~s (* = Fl [s] *) 



(* part B. integrand has a pole at s = (2 + each zeta zero) *) 
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(* set B[s] = Zeta[s - 2]; set A[s] = everything else *) 

A[s] = Zeta[2 s] Zeta[2 s - 2] Zeta[2 s - 4] / (Zeta[s] Zeta[s - 1] ~2) * x"s/s ; 
B[s] = Zeta[s - 2] ; 

residueFormula[l , s]/x"s /. s -> s+2 (* = F2 [s] *) 

(* part C. integrand has a pole of order 2 at s = (1 + each zeta zero) *) 
(* set B[s] = Zeta[s - 1]"2; set A[s] = everything else *) 

A[s] = Zeta [2 s] Zeta [2 s - 2] Zeta [2 s - 4] / (Zeta[s] Zeta[s - 2] ) * x's/s ; 

B[s] = Zeta[s - 1] ~2 ; 

expr = residueFormula[2, s] ; 

(* we will evaluate expr only at s - 1 = (zeta zero), so remove all Zeta[s-1] *) 
expr2 = Together [expr /. Zeta[s-1] -> 0]; 

(* compute the residue at s = (1 + each zeta zero), so make this substitution *) 
exprS = expr2 / . s -> s+1 ; 

numer = Collect [Numerator [exprS] , Log[x] , Simplify]; 
numer[[2]]/x-(s+l) (* = F3[s] *) 

numer [[l]]/(x" (s+1) Log[x]) (* = F4[s] *) 
Denominator [expr 3] (* = G[s] *) 

Note the presence of C(2s) in Fi, F2, F3, and F4. If s = —2k is a real zeta zero, then so is 
2s = —Ak. Therefore, the sum over the real zeta zeros vanishes. 



N = 50 N = 100 




Figure 48. Approximating the sum of X{n)a{ny using equation (22.2) with 
iV = 50 and iV = 100 



Figure 48 shows the summatory function of X{n)a{n)'^ and the approximations that come from 
equation (22.2) using the first 50 pairs and 100 of zeta zeros. Although both approximations 
can reproduce the large swings in the step function, neither approximation is particularly good 
at reproducing the details of the step function. Further, the approximation does not appear to 
be improved by using 100 pairs of zeros instead of 50. (As usual, using more zeros does enable 
the sum to reproduce the general shape of the step function for larger x). On the other hand, 
Figure 49 shows that the integral approximation (22.1) works quite well. 
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c = 3.5,T= 144 



c = 3.5, T = 237 



10000 



5000 



-5000 



-10000 




10000 - 



-5000 - 




Figure 49. Approximating the sum of X{n)a{ny using the integral in (22.1) 

23. Tallying X{n)a{n)T{n) 
This Dirichlet series for X{n)a{n)T{n) holds for s > 2 [9, eq. D-46], [20, p. 232] 

A(n)ff(n)T(n) C(2s)'C(2s - 2)^ 



E 

n=l 



The summatory function is 



c(2s - msms - ir 



T{x) = \{n)(7{n)T{n). 



k=l 



The integral approximation to the modified summatory function To(x) is 



c+iT 



C(2s)2C(2s-2)2 



ds. 



The sum over the residues, which we hope will approximate the summatory function, is 



(23.1) 



To{x) ~ oia/x + a2Va; log(a;) + asx^^'^ + 040;^/^ log( 



X) 



(23.2) 



.Fi{pk) + F5 (pfc) log(a;) 



G2 (Pk) 



where 



C i-l) (C il) (-241og(A) + 3 + log(27r)) + C (i)) + C H) C H) - 27C (D C (-^) 



where 7 is Euler's constant and A ~ 1.282427 is Glaisher's constant defined by 

1 

~ 12 
1 



0.321773, 



log(A) 



C'(-i), 



02 



-0.075348, 
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2C(3) 
3-^(0 C(i) 



-3c(Dc(3)c'Q-cQ(3C(3)c'(^).c(^)(c(3)-oa^^^ 



~ -0.002403, 

04 = ^ ~ 0.010059, 

F2(s) = -C(2s)C(2s-2) 

- l)C(2s)C(2s - 2)C(2s - l)C"(s) + C'(s) 

(2<(25)C(25 - 2)C(2s - 1)C'(^ - 1) + as - 1) 
(C(2s) (C(2s - 2) (C(2s - 1) + 2<'(2s - 1)) 

- 4<(2s - l)C'(2s - 2)) - 4<(2s - 2)C(2s - l)C'(2s) 

F3(5) = <(5 - l)C(2s)\(2s - 2)^(25 - l)C{s), 

F4(5) = C(2s)C(2(s + l)) 

( - + l)C(2s)C(« + l)C(2(s + l))C(2s + i)C"(^) - C{s) 
(C(2s) (2{s + l)C(2(s + l))C(2s + l)C{s + 1) + C(s + 1) 

(C(2(8 + 1)) (C(2s + 1) + 2{s + l)C{2s + 1)) - 4(s + l)C(2s + l)C{2{s + 1)))) 
- 4(s + l)C(s + l)C(2(s + l))C(2s + l)C'(2s) 



i^5(«) ^{s+ m^sms + i)c(2(5 + i)n{2s + i)as), 
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and, finally, 

The following Mathematica code will perform these calculations: 

(* part A. integrand has a pole at each zero of Zeta[2s - 1] *) 
(* set B[s] = Zeta[2s - 1]; set A[s] = everything else in the integrand *) 
A[s] = Zeta[2s]"2 Zeta[2s - 2] ~2 / (Zeta[s]''2 Zeta[s - 1] "2) * x"s/s ; 
B[s] = Zeta[2s - 1] ; 

(* compute the residue at s = (1 + (zeta zero))/2, so make this substitution *) 
residueFormula[l, s] / x"s /. s -> (1 + s)/2 (* = Fl [s] *) 

(* part B. integrand has a pole of order 2 at each zero of Zeta[s] *) 
(* set B[s] = Zeta [s] "2; set A[s] = everything else in the integrand *) 
A[s] = Zeta[2s]''2 Zeta[2s - 2] "2 / (Zeta[2s - 1] Zeta[s - 1] "2) * x~s/s ; 
B[s] = Zeta[s] "2 ; 
exprB = residueFormula[2, s] ; 

(* we will evaluate exprB only at s = zero of zeta function, so remove Zeta[s] *) 

exprB2 = Together [exprB /. Zeta[s] -> 0]; 

numerB = Collect [Numerator [exprB2] , Log[x] , Simplify]; 

nuinerB[[2]]/x"s (* = F2 [s] *) 

numerB [[l]]/(x"s Log[x]) (* = F3 [s] *) 

Denominator [exprB2] (* = Gl [s] *) 

(* part C. integrand has a pole of order 2 at each zero of Zeta[s - 1] *) 
(* set B[s] = Zeta[s-1]"2; set A[s] = everything else in the integrand *) 
A[s] = Zeta[2s]~2 Zeta[2s - 2] ~2 / (Zeta[2s - 1] Zeta[s]"2) * x"s/s ; 
B[s] = Zeta[s - 1] ~2 ; 
exprC = residueFormula[2, s] ; 

(* we evaluate exprC only at s-1 = zero of zeta function, so remove Zeta[s-1] *) 
exprC2 = Together [exprC /. Zeta [s-1] -> 0]; 

(* we want the residue at s = (1 + each zeta zero) ; so make this substitution *) 
exprCS = exprC2 /. s -> s+1; 

numerC = Collect [Numerator [exprCS] , Log[x] , Simplify]; 
numerC[[2]]/x"(s+l) (* = F4[s] *) 

numerC[[l]]/(x"(s+l) Log[x]) (* = F5[s] *) 
Denominator [exprCS] (* = G2[s] *) 

What about the sum over the real zeta zeros? This sum would have the same form as the 
terms in (23.2), but with replaced with —2k. But because of the C(2s) that is present in F2, 
F3, F4, and F5, we see that those parts of this sum would be zero. However, the part of the 
sum involving Fi{s) evaluated at s = —2k, that is, 

AI 

x^^Fi{-2k) 

k=l 

is quite badly behaved: the coefficients Fi{—2k) rapidly grow large. The first ten are ap- 
proximately 0.56426, -0.96909, 6.478, -90.794, 2197.1, -82446, 4.4633 ■ 10^ -3.3062 ■ 10^ 
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3.2185 ■ 10^°, and —3.9889 ■ 10^^. As we saw in section 18, if the sum over real zeros does not 
improve the accuracy of the approximation to T{x), we may be able to omit those terms. All 
we have to do is integrate around a rectangle in the complex plane that does not enclose the 
real zeta zeros. 



N = 50 



N= 100 



1500 - 



1000 - 



500 - 



-500 - 



-1000 - 




1500 - 



500 - 



-1000 - 




Figure 50. Approximating the sum of X{n)a{n)T{n) using (23.2) with = 50 
and A^ = 100 




Figure 51. Approximating the sum of X{n)a{n)T{n) using the integral in (23.1) 



Figure 50 shows the summatory function of X{n)a{n)T{n) and the approximations that come 
from equation (23.2) using the first 50 pairs and 100 of zeta zeros. Neither sum is particu- 
larly good at reproducing the details of the step function. Figure 51 shows that the integral 
approximation (23.1) works quite well. 

24. Tallying Greatest Common Divisors 
Define the function P{n) to be 

n 

P{n) = gcd(fc, n). 

k=l 
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The following Dirichlet series for P{n) holds for s > 2 [35, eq. 15]: 

" P(n) C(-5-l)" 



E 



C(s) 



n=l 

This Dirichlet series comes from the fact that P{n) is the convolution of {(f) * Id){n) where 
Id{n) = n is the identity function, so the dirichlet series for P is the product of the Dirichlet 
series of (p and Id: 

^Pjn) _ as-l) ^n _ as-i) _ 

71=1 ^ \ / j^ — l -3 \ J 

Here, we used the Dirichlet series for 0, given by (11-1). 
The summatory function is 



T(x) = ^P( 



[n) 

k=l 

The integral approximation to the summatory function T{x) is 

Ux)^^. ^-^-ds. (24.1) 

The integrand has a pole of order 1 at s = 0, and a pole of order 2 at s = 2. There is also a 
pole at each zeta zero. The residue at s = is —1/72. The residue at s = 2 is 

3x^ ((47 - l)7r2 - 12C'(2)) Sx^ log(a;) 



+ 



27r4 7r2 
The sum of the residues, which we hope will be close to To(x), is 

N ^ , . N 2 ^ 



To(x) ^ a, + a.x^ + a.x^ log(x) + 25? ( ^^^^^ ) + ---^^y^;^ (24.2) 

where 



PkC'{Pk) / (-2A;)C'(-2A;) 



ai = -1/72 ~ -0.013889, 

3((47- l)7r2- 12C'(2)) 

02 = ~ 0.37217, 

27r 

and 

as = 3/7r2 ~ 0.30396. 

Figure 52 shows the graph of the sum approximation (24.2) without using any zeta zeros. 
Figure 53 shows the graphs of the sum approximation (24.2) using 50 and 100 pairs of complex 
zeta zeros. Figure 54 shows the integral approximation (24.1). 

The sums using 50 zeta zeros is pretty poor. The obvious thing to try is to use more zeros. 
But that only makes the graph worse. If one makes a table of values of, say, the first 50 
coefficients in the first sum in equation (24.2), that is, 

C(P.-I)^ 
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Figure 52. Approximating the sum of P{n) with the first three terms of (24.2) 
N = 50, M = 



N = 100, M = 




Figure 53. Approximating the sum of P{n) using 50 and 100 pairs of zeta zeros 

for 1 < A; < 50, one sees that neither the real nor the imaginary parts of the coefficients are 
getting smalL This is in contrast to previous cases where the sum yielded a good approximation. 
It may be that the first sum in (24.2) does not even converge. 
To see what's going on, let's compare the two functions 

gpk - 1) 

pkC'ipk) 

from (11.3), the sum for the summatory function of (f){n), and 

PkC'ipk) 

from the summatory function of this section. In particular, let's graph the function 

C(^ - 1) 
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c = 2.5, T = 144 




Figure 54. Approximating the sum of P{n) using an integral 




Figure 55. Left: function (24.3); right: function (24.4) 



where s lies on the hne 1/2 + it, that is, 

C(l/2 + it- 1) 
(1/2 + it)C{l/2 + it) 

for, say, < t < 100. Let's also graph 

C(l/2 + ^t-l)^ 
(1/2 + it)C{l/2 + it) 

over the same interval. 



(24.3) 



(24.4) 



It is also stated [5, Theorem 1.1] that, for x > 1 and every e > 0, 

T{x) = a2X^ + a^x^ log(x) + 0,(x^+^+') 

where 02 and 03 have the values given earlier in this section, and where 1/4 < ^ < 131/416 is 
the exponent that appears in the Dirichlet divisor problem. 
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25. Powerful ( "power- full" ) Numbers 

We now consider integers that, in some sense, are the opposite of squarefree integers. Let k 
be an integer > 2. If the prime factorization of n is 

then we say that n is fc-full if every exponent aj > k. (We also define 1 to be fc-fuU for all 
k > 2). The characteristic function of the fc-fuU numbers is 



fkin) 



1 if n is fc-full, 
otherwise. 



If = 2, then we call the 2-full numbers "square-full", and we have this Dirichlet series [15, 
p. 33]: 



Mn) _ C(2s)C(35) 

n=l ^ ' 



E 



which converges for s > 1/2. 

The summatory function, which counts the square-full numbers up to x, is 



n=l 

The integral approximation to the summatory function Tq{x) is 

The integrand has poles of order 1 at s = 0, at s = 1/3, and at s = 1/2, where the residues 
are -1/2, 

C(2)' 

and 

(i) 
C(3)' 

respectively. 

From these, from the residues at each zero of C(6s), that is, at s = p/6, we get the following 
sum, which we hope approximates the summatory function Tq{x): 



where ai = —1/2, 



as = ^ ~ -1.48795, 
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and 



as 



c(l) 

C(3) 



2.17325. 



N = 0, M = 




Figure 56. Counting squarefuU numbers using (25.2) with N = M = Q 



N = 50, M = 3 



c = 1,T= 144 





Figure 57. Counting squarefull numbers using (25.2) (left) and (25.1) (right) 



Figure 56 shows the summatory function and the approximation that comes from setting 

1/2. The left side 



1/3 



N = M = Q m. equation (25.2), that is, using only the terms Oi + 
of Figure 57 shows the approximation using = 50 pairs of complex zeta zeros and M = 3 real 
zeros. Without including the 3 real zeros, the approximation would be visibly, but slightly, too 
low. The right side of Figure 57 shows the approximation based on the integral approximation 
(25.1). Figure 58 shows the same approximations, but for a larger range of x. 
Bateman and Grosswald [4] proved that the number of squarefull numbers up to x is 

a/3C(2/3) 



X 



C(2) 



C(3) 
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200 300 400 500 



100 200 300 



Figure 58. Counting squarefull numbers for x < 500 

where a is a positive constant, and 

u{x) = log(x)^/^log(log(x))-3/7. 

The approximation in equation (25.2) with N = M = consists of the first two terms of 
this estimate, plus the constant —1/2, but we have not proved that the difference between our 
estimate and the summatory function is smalL 



26. What's Really Going on Here? An Analogy With Fourier Series 

When you stop to think about it, Fourier series are pretty amazing. By taking a hnear 
combination of a countable number of "wavy" functions (sines and/or cosines), we can get a 
function that is a straight line over an entire interval that contains an uncountable number of 
points! 

Fourier series do have one limitation, however: they are periodic. That is, the sum of sines 
and cosines will represent a desired function over only a finite interval, say, between and tt, 
at which point the sum will repeat. 

We will see that, like Fourier series, the sums we computed above are also linear combinations 
of "wavy" functions. But in one sense, our sums are even more amazing than Fourier series. 
The step functions we are trying to approximate are quite irregular and are not periodic. 
Nevertheless, our sums often appear to approximate the step functions all the way out to 
infinity! 

Consider the sums we've obtained over pairs of complex zeta zeros, for example, (8.6), (9.3), 
(10.2), (11.3), (12.2), or (13.3). 

A typical term in one of these sums has the form Ckx'''' where Ck is a complex coefficient, say 
Ck = CLk + ibk- Let's write the fc*'* complex zeta zero as pk = 1/2 + itk. Then 



3,1/2 . ^itk ^ ^1/2 . exp(itfc log(a;)) 



X 



1/2 



[cos{tk log(x)) + i sin(tfc log(x))]. 



70 ROBERT BAILLIE 

Then 

Cfcx''* = x^/^ ■ (ak + ibk) ■ [cos(tfc log(x)) + ism{tkhg{x))]. 
We will want the real part of this product, which is: 

3f?(cfcx'''=) = x^/^ ■ [afeCos(tfclog(x)) - 6fc sin(tfc log(a;))] 

= ttk ■ x^l'^ cos(tfc log(x)) - hk ■ x^l'^ sin(tfc log(a;)). 

To make this more concrete, let's consider the first such term, with pi ~ 1/2 + 14.135i, and, 
for comparison, the tenth such term, with pio ~ 1/2 + 49.774i, and let's take the coefficients C\ 
and Cio to be 1. Then the real parts of the products are 

x^/^ cos(14.135 log(x)) 

and 

x^l^ cos(49.7741og(x)). 

Here are the graphs of these two functions. Both curves are bounded by the envelopes ±x^/^. 




Figure 59. x^/2cos(14.1351og(x)) and x^/^ cos(49.7741og(x)) 



Again, the approximating sums are linear combinations of functions that look like these. We 
can think of the functions 

x^l'^ cos(tfc log(x)) 

and 

x^*^^ sin(tfc log(x)) 

as hasis functions for our approximating sums, in the same sense that sin(nx) and cos(nx) are 
the basis functions for Fourier series. 

This analogy also helps explain why there appears to be a Gibbs-type phenomenon in some 
of the graphs above. 
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27. Counting Primes in Arithmetic Progressions 

In section 4, we saw how to use the first pairs of zeros of the Riemann zeta function to 
count the primes up to a; . Here, we will make a first pass at using zeros of Dirichlet L-functions 
to count primes in arithmetic progressions. 

Let q and a be integers with g > 0. If g and a are relatively prime, then the arithmetic 
progression qn + a where n = 0, 1, 2, 3, . . . contains an infinite number of primes. This is 
Dirichlet's Theorem [18, p. 76]. The number of such primes that are less than or equal to x is 
a function of q, a, and x that is usually denoted by nq^a{x). For an introduction to characters 
and Dirichlet L-functions, see [18, Chapter 6]. 

Equations (27.1) through (27.4) below come from equations (6) through (8) in [19]. See [29] 
for an estimate of the error term for these equations. 

a/2 



= TTIT +^(^'^'«)77ITT7Z7ZT (27-1) 



-— + E{x,q,a)-—— 
(j){q) (f){q) log(x) 

where 

Fix, q, a) = -c{q, a) - xiayE{x, x), (27.2) 

Xt^XO 

c(g, a) = —1 + ^(square roots of a, mod q), (27.3) 

and 

E(x,x) = Y ^ • (27.4) 



k=l 



In equation (27.1), the main term is Except for c{q,a), the rest of equation (27.1) 

depends on zeros of L-functions. 

In equation (27.2), the asterisk denotes the complex conjugate, and we sum over all characters 
X (mod q), except for the principal character, xo- 

In equation (27.3): if a has no square roots (mod q), that is, if a is a quadratic nonresidue 
(mod g), then c(g, a) = —1. Otherwise, a must have at least two square roots (mod q), which 
makes c(g, a) > 1. This causes a bias in "prime number races": Hq^aix) often exceeds 7rq^b{x) if 
6 is a square (mod q) and a is not. For example, 773^2(2;) > vr3^i(x) for all x < 608981813029. 
See [11] for more details on this fascinating topic. 

In equation (27.4), for the character x (mod q), we sum over complex zeros of the corre- 
sponding Dirichlet L-function. We include 2N such zeros in the sum: namely, the N zeros 
having the smallest positive imaginary part, and the N having the least negative imaginary 
part. Note that these zeros do not necessarily occur in conjugate pairs, pk is the imaginary 
part of the k^^ zero of the L-function that corresponds to the character x- 

The solid curve in figure 60 is 7r3^i(x). The dashed curve is an approximation to 773^1 (x) using 
equations (27.1) - (27.4) with no Dirichlet L zeros taken into account. That is, the dashed 
curve is the approximation 
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Figure 60. 7r3 i(a;) (solid), and an approximation (27.5) using no zeros of L- 
functions (dashed) 



7r(x) c{q,a) x^l"^ 



0(g) 0(g) log(a;) 

(This approximation looks like a step function in spite of the second term on the right, 
because the second term is relatively small and doesn't vary much for x between 2 and 50). 
Note that this approximation uses vr(a;), so the dashed curve jumps up at every prime. 




Figure 61. 7r3^i(x) and an approximation using 200 zeros of L-functions 

Figure 61 shows 7r3^i(x) and an approximation using the first 2A^ = 200 zeros of L-functions 
in equations (27.1) - (27.4). Note the prominent glitches at x = 17, 23, 29, 41, and 47. These 
are places where 7r(a;) jumps, but 7r3^i(x) does not. 

Warning: there is a bit of empirical fakery in this graph! The pair of arithmetic progressions 
{3n + 1, 3n + 2} omits one prime, namely, 3. The value for 7r3 i(x) that comes out of equation 
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(27.1) is too large by about 1/2. The same is true for 7r3^2{x)- So, to account for the missing 
prime among the two progressions {3n + l,3n + 2}, the values shown in figure 61 are 1/2 
less than the values we get from equation (27.1). In spite of this legerdemain, the sums over 
zeros of L-functions provide a decent approximation to the counts of primes in these arithmetic 
progressions. 




Figure 62. vr3_2(x) and an approximation using 200 zeros of L-functions 

Figure 63 shows 7r4_i(x) and 7r4^3(x) and their approximations using 2N = 200 zeros of L- 
functions. In order to account for the fact that this pair of progressions omits the prime 2, the 
values shown in the graphs are 1/2 less than the values we get from equation (27.1). 




Figure 63. 7r4,i(x) (left) and vr4,3(x) (right) 



In figures 64 and 65, we use 2N = 200 zeros of L-functions to approximate 7rio,i(x), 7rio,3(x), 
7rio,7(x), and 7rio,9(x), the four progressions mod 10, up to x = 200. These four progressions, 
taken together, omit two primes, namely, 2 and 5. For the progressions with a = 1, 3, 7, 9, we 
adjust the values by, respectively, —1/4, 0, —3/4, and —1. These are empirical guesses. It is 
not clear why these values seem to work better than, say, —1/2 for all of the progressions. 
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Figure 64. 7rio,i(x) (left) and 7rio,3(x) (right) 




Figure 65. 7rio,7(x) (left) and 7rio,9(x) (right) 



28. Questions For Further Research 
The answers to the following questions are not known to this author. 

Can the Riemann-von Mangoldt formula (equations (4.1) - (4.2)), which counts primes using 
a sum over zeta zeros, be obtained by applying Perron's formula? 

In section 15.2 that considers the sum of r(n^), when we use 200 zeros in the sum, the 
graph seems worse in some ways than when we use only 100 zeros. The same thing happens 
in section 16 with a{n)T{n), in section 23 with X{n)a{n)T{n), and in section 24 where we tally 
CCD's. 

Why does this happen? Can we predict in which cases will this occur? Would the approxi- 
mation get worse and worse if we used even more zeros? 

In section 18 on A(n)2'^*^"\ section 19 on A(n)r(n^) and section 20 on A(ra)r(n)^, and 23 on 
X{n)a{n)T{n), if we extend the rectangle to the left to include some real zeros of zeta, the sum 
of the residues appears to diverge. Can we prove this? For which Dirichlet series will this 
happen? 

In section 26 on the analogy with Fourier series: is there anything more we can say that's 
useful or interesting? 



EXPERIMENTS WITH ZETA ZEROS AND PERRON'S FORMULA 



75 



In section 27 on counting primes in arithmetic progressions, why, for q = 10, do we need 
different offsets for a = 1, 3, 7, and 9? 

In several cases, there are provable theorems which say (more or less) that a sum over zeta 
zeros does, in fact, approach the respective arithmetical function. These cases are: 

• equations (4.1) - (4.2) (counting primes), 

• equation (8.6) (computing ip{x)), 

• equation (9.5) (computing the Mertens function), 

• equation (10.3) (counting squarefree integers) (a similar result holds for cubefree inte- 
gers), 

• equation (13.4) (counting squarefree divisors). 

Can similar theorems be proven for the other arithmetical functions that are considered here? 
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